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ANALYSIS  AND  EVALUATION  OF  SUPPRESSIVE  SHIELDS 


I.  INTRODUCTION 

Tliis  final  technical  report  documents  work  performed  under  Contract  DAAAl  5-75-C-0083 
for  Edgewood  Arsenal  in  support  of  the  suppressive  structures  program.  Principally,  it  includes 
work  performed  during  the  time  period  of  January  1976  through  June  1977;  however,  some 
earlier  work  has  been  included  in  summary  form.  Also,  some  work  performed  during  this 
time  period  and  documented  in  separate  letter  reports  has  not  been  included.  Work  performed 
prior  to  January  1976  is  documented  in  References  I through  9. 

Contents  of  this  report  cover  three  different  aspects  of  the  work: 

• The  use  of  energy  methods  to  predict  deformations  in  ulast-loaded  structures, 

• Strain  data  analysis  for  the  Category  1 1 /4-scale  model  tests,  and 

• Calculation  of  pressure-time  histories  produced  by  burning  propellant  in 
vented  enclosures. 

Chapter  II  contains  recent  work  on  the  development  of  energy  solutions  for  structural 
response.  This  recent  work  includes  a study  of  the  influence  of  hie  deformed  shape  on 
the  accuracy  of  the  solution  obtained,  energy  solutions  for  coupled  response,  energy  solu- 
tions for  combined  elastic-plastic  behavior  of  beams  and  strings  and  the  construction  of 
general  graphical  solutions  for  blast-loaded  beams.  Chapter  III  contains  a summary  of 
solutions  which  have  been  developed  over  the  total  contract  period. 

Analysis  of  strain  data  from  the  Category  1 1 /4-scale  model  tests  is  covered  in  Chapter 
IV.  Comparisons  are  made  between  measured  strains  and  strains  predicted  by  the  approxi- 
mate energy  methods  and  by  finite-element  methods.  Chapter  V describes  the  calculation 
procedures  and  gives  results  for  M 10  propellant  burning  in  vented  enclosures. 


II.  ADDITIONAL  SOLUTIONS  WITH  ENERGY  METHODS 
A.  Importance  of  the  Assumed  Deformed  Shape 

In  using  energy  solutions  to  compute  maximum  deformations  or  strains  in  blast  loaded 
structural  components,  we  select  a deformed  shape  with  appropriate  boundary  conditions. 
Usually,  this  assumed  deformed  shape  is  either  the  first  mode  from  an  infinite  series  of  modes, 
or  the  static  deformed  shape  carried  over  to  a dy.namic  analysis.  We  will  demonstrate  that 
either  assumption  can  give  excellent  predictions  of  strains  or  deflections  using  eiV/ier  assumed 
deformed  shape.  This  generalization  pertains  for  both  elastic  and  plastic  response  of  struc- 
tural components,  and  is  true  for  both  quasi-static  and  impulsive  transient  loads. 
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I.  Influence  of  Higher  Modes 


A general  deformed  shape  for  a simply-supported  elastic  beam  loaded  with  a 
uniform  quasi-static  pulse  is  given  by  Eq.  (I ). 


' ~ sin 


Nnx 


(1) 


N=  1,3.5 


The  even  modes  are  missing  because  of  symmetry,  and  cosine  contributions  are  missing  be- 
cause of  simply-supported  boundary  conditions.  For  bending  only,  the  strain  in  the  beam  is 


A'^y  v^woh  ^ 

^ ~ ~ 9.^  ^ 

^ N=1,3,S 


sin 


N-nx 


In  a rectangular  beam,  t 


'e  is  given  by; 

W/2  «/7, 

Vol  =4b  Jf  Ah  j Ax 


(2) 


(3) 


The  strain  energy  U is: 


f f dK, 

y ^ 


(4a) 


or: 


f/  = 


Iv*  Ewlb 


g4 


21.  /M  \ 

--  j h'^Ah  S 


N=1.3,5 


(4b) 


The  cross-products  associated  with  squaring  the  strain  integrate  to  zero  in  the  preceding 
equation  because  of  orthogonality  of  the  mode  shapes.  Performing  the  required  double 
integration  gives  for  the  strain  energy; 


^ ;vMl 


48  6=* 


(5) 


N=l.3,5 


Next  we  compute  the  maximum  possible  work  imparted  to  the  structure  with  the  deformed 
shape  of  Eq.  ( 1 ).  This  work  equals; 
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oo 


8/2 

wk-2pbw„  S Afj  f 
N=l.3,5  0 


. NitX 

Sin d-v 

il 


Or,  after  integrating: 


, 2pfcw»„  C 
wk  = — ■ — - 

TT 


E 

N=  1 ,3,5 


Ajv 

N 


(6) 


(7) 


The  Rayleigh-Ritz  method  can  be  used  to  obtain  the  relative  amplitudes  of  the  different 
modes.  This  is  accomplished  by  subtracting  the  work  from  the  strain  energy, 

OO  OO 


( U - wk)  = 


48  gs 


N^Al  - 


IpbWpfi 


A/^ 

N 


N = 1 ,3.5 


N=  1 ,3,5 


(8) 


and  differentiating  with  respect  to  Afj  so  the  energy  difference  is  minimized.  This  procedure 
yields: 


d{U-wk)  n^EwlbH^  . IpbyVo'i  1 
3/1  A-  " 24  2*  ^ It  N~ 


or: 


As  - 


A'* 


(10) 


If  one  substitutes  As  from  Eq.  (10)  back  into  the  preceding  equations,  a series 
solution  for  the  deformed  shape  is  obtained.  The  influence  of  higher  modes  can  be  obtained 
by  dividing /1a  by  the  amplitude  of  the  first  mode.  This  step  yields 


As  _ I 

A^~ 


(11) 


For  a three-mode  solution,  the  deformed  snape  will  equal: 


>'  = mwo 


■nx 

sill  — + 
2 


1 

243 


3itx  1 57r.r 

sin 1 — sin 

2 3125  2 


(12) 


The  parameter  iii  equals  1 .0038  if  y is  to  equal  at  mid-span  (at  x/S.  to  0.5).  Substituting 
for  tn,  differentiating  Eq,  (12)  twice.,  and  multiplying  by  —fi  gives  the  strain. 


' . Ttx  1 37rA  1 5nx 
siUYT-sin  — -f—sin- 


(13) 
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From  Eq.  (4a), 


1.0076 


0 


Performing  the  double  integration  gives  for  the  strain  energy: 


U=  1.053 


(v*EbwlH^\ 
\ 48  8*  / 


(15) 


The  factor  1 .053  has  been  separated  from  the  other  terms  in  the  precedir;  -quation, 
as  this  5.3  percent  increase  in  the  strain  energy  U is  the  effect  of  adding  the  third  and  fifth 
mode  contributions  to  the  first  mode  estimate.  A similar  procedure  will  be  followed  in 
presenting  the  influence  of  high  modes  in  subsequent  relationships. 


The  work  which  must  be  equated  to  this  strain  energy  is  given  by: 


8/2 

wk  ~ 2pb{  1 .0038^0 ) J 
0 


TTX 

sin  — + 
t 


1 

243 


. 37r;c 

sin 

2 


+ 


I . SttxI 
3125  2 J 


(Lr 


(16) 


or,  after  integrating: 


1.00524 


(17) 


The  influence  of  higher  modes  on  the  work  is  less  than  a 1 percent  increase. 
Equating  the  work  to  the  strain  energy  and  solving  for  a nondimensionalized  mid-span 
deformation,  w„///8^,  yields: 


WqH 

8» 


0.9546 


96p8* ' 


(18) 


This  equation  implies  that  the  actual  deformation  will  be  less  than  that  estimated  with  the 
one  mode  approximation,  but  the  difference  is  less  than  5 percent.  The  maximum  strains 
can  also  be  estimated.  Substituting  H/2  for  h,  Eq.  ( 1 8)  into  Eq.  ( 1 3),  and  2/2  for  x gives 
for  the  maximum  strain: 


0.9304  (19) 

Equation  (19)  shows  that  the  one  mode  approximation  also  overesEmates  the  strain  in  this 
illustration,  but  that  this  error  is  le'^s  than  7 percent. 
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We  present  one  more  illustration  to  show  that  these  same  conclusions  can  be 
reached  for  plastic  as  well  as  elastic  response,  and  for  extensional  as  well  as  bending  behavior. 
As  the  second  illustration,  consider  a string*  with  the  deformed  shape  given  by  Eq.  (1).  ITie 
string  has  a different  shape  from  a beam;  ho'vever,  this  difference  is  reflected  in  the  An 
coefficients.  In  a string, 


_ /dy  \ * 
' 2 \d)c/ 


^ = ^.1:  /VM^cos^-—  (20 

N=  1 .3.5 

S 

The  plastic  strain  energy  equals  the  yield  force  times  the  strain  integrated  over  the  length: 


N=  1 .3.5 


/ cos^  (^)  ^ 


or,  after  integrating: 


42 


N-l.3.5 


If  the  loading  is  a quasi-static  pressure  of  intensity  p,  then  the  work  is  given  by  the 
previously  obtained  Eq.  (7).  We  again  use  tne  Rayleigh-Ritz  method  to  obtain  amplitudes  of 
the  different  modes.  Although  our  system  is  not  conservative,  it  is  linear  because  we  have 
chosen  to  consider  rigid-perfectly  plastic  mateiial  behavior  over  the  full  range  of  deformation. 
Proceeding  as  before  by  taking  the  difference  in  the  strain  energy  and  work,  we  have 


{U  — wk)  = 


ir^  Fy  wS 


N^Al  - 


Ipbwp  £ 


N = 1 ,3.5 


N=l,3.5 


Differentiating  with  respect  to  An  and  setting  the  result  equal  to  zero  yields  ior  An  - 


= (^P^\  _L 


*A  siring  is  defined  here  as  an  element  with  negligible  bending  stiffness. 
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Hence,  the  amplitude  of  various  modes  relative  to  the  first  one  is 

(25) 

/4 , N* 

If  we  now  proceed  with  a dynamic  solution,  the  deformed  shape  for  a three  mode 
solution  is  given  by: 


y = mv/o 


r tta:  1 


3itx  1 

sinT  + T^^'n  — 


The  parameter  m equals  1 .0299  because  y is  equal  to  Wq  at  mid-span.  Substituting  into 
Eq.  (20a)  gives  as  an  approximation  for  the  strain: 


e = 1 .0607 


The  plastic  strain  energy  is 


.28*  J L 


nx  1 3»rjc  1 5^x1  * 

cos  — + — cos  -j-  + — cos  — 
£9  £25  2 


U=2  f Fy€dx 


U=  1.0607 


[^]  / [< 


irx  1 3rx  1 5vr.xl*  , 
cos  — + — cos  f-  “ cos  — ' djr 
£ 9 £ 25  £ J 


Performing  the  desired  integration  yields: 


L 4£  J 


Tlie  work  which  must  be  equated  to  this  strain  energy  is  given  by: 


e/2 

t-  [ tx  \ , 3irx  1 . Svx~ 

wk  = 2pi»(1.0299wo)  J sin  Y + ^ "T"  Ti?  "T 

0 ^ 


or,  after  integrating, 


wk=  1.0443 
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Equating  the  work  to  the  strain  energy  and  solving  for  a nondimensionalized  mid-span  defor- 
mation. Wo/fi,  gives; 


e 


0.8553 


8pf>g 

■n^Fy 


(32) 


This  plastic  extensional  solution  is  not  as  exact  using  only  one  mode  as  the  elastic 
bending  solution.  Nevertheless,  the  deformation  is  predicted  to  within  1 5 percent,  and  the 
single  mode  does  overestimate  the  deformation.  The  maximum  strains  can  also  be  compared. 
This  maximum  occurs  at  .v  = 0.  Substituting  Eq.  (32)  into  Eq.  (27)  and  .settings  - 0 gives 
for  the  maximum  strains: 


tw  = ■1.0282 


(33) 


As  is  seen  in  Eq.  (33),  the  strain  is  predicted  fairly  accurately  using  a one  mode 
approximation.  The  difference  between  a one  mode  (the  first)  and  a three  mode  approxi- 
mation is  less  than  3 percent.  Considering  the  mathematical  complexities  which  are  intro- 
duced by  adding  more  modes,  the  additional  effort  is  probably  not  worth  it. 


2.  Influence  of  Other  Shapes 


To  illustrate  the  influence  of  assumed  deformed  shape  on  results,  we  analyzed 
elastic  bending  and  plastic  bending  of  simply-supported  beams  loaded  with  a uniform  impulse. 
Examples  of  these  solutions  are  given  in  Part  C of  this  section  and  will  not  be  repeated  here. 
Equations  are  summarized  in  Chapter  III. 

For  an  clastic  analysis,  we  evaluate  computed  results  for  three  different  deformed 
shapes.  The  first  shape  is  a parabola,  the  second  is  the  first  mode  sine  wave,  and  the  third 
shape  is  the  static  one  for  a uniformly  applied  load.  Because  each  of  these  deformation 
patterns  must  satisfy  the  essential  boundary  conditions  of  a simply-supported  beam,  they 
will  be  similar  in  shape.  Nevertheless  they  represent  three  distinct  deformation  patterns. 

All  three  solutions  give  similar  results  for  the  strain  energy  17,  the  maximum  defor- 
mation Wo , and  the  maximum  strain  Cmax  i as  can  be  obwrved  in  Table  I.  The  results  in  Table 
I have  been  nondimensionalized  so  numerical  coefficients  can  be  compared  directly.  Only  the 
numerical  coefficient  differs  slightly  in  all  of  these  solutions.  Both  the  first  mode  and  static 
deformed  shape  approximations  yield  the  same  answer  to  within  a few  percent. 

The  second  group  of  illustrations  includes  plastic  bending  in  a beam  loaded 
dynamically  with  a uniform  impulse.  One  more  deformed  shape  will  be  added  to  the  results 
for  the  three  elastic  shapes.  The  fourth  shape  is  a static  hinge  yielding  plastically  in  the  center 
of  the  beam  while  the  rest  of  the  beam  remains  rigid.  This  deformed  shape  is  a common  one 
in  civil  engineering  applications,  Whereas  the  other  three  deformed  shapes  distribute  the 
deformation,  the  stationary'  hinge  concentrates  the  deformation.  Table  II  compares  non- 
diniensionalized  numerical  coefficients  for  strain  energy,  maximum  plastic  deformation, 
and  maximum  plastic  strain  for  these  four  different  deformed  shapes. 
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TABLE  I IMPULSIVE  BENDING  SOLUTION  FOR  A SIMPLY-SUl'i'OKTEO  ELASTIC  BEAM 


In  the  plastic  beam,  the  first  mode  and  static  defomried  shape  give  deformations  and 
strains  which  are  close  to  one  another.  The  stationary  hinge  yields  no  rational  estimate  for 
strain  because  there  is  no  gauge  length  associated  with  a concentrated  hinge.  The  deformation 
is  also  much  larger  for  a stationary  hinge  than  comparable  deformations  from  distributed 
deformed  shapes.  A concentrated  hinge  yields  much  too  small  a strain  energy  and  too  large 
a deformation.  Distributed  deformed  shapes  are  much  closer  to  reality  and  give  more  accurate 
predictions. 
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As  tlK‘so  illustrations  show,  either  a first  mode  approximation  or  a static  deformed 
shape  is  a good  approximation.  We  would  recommend  a first  mode  approximation  for  sym- 
metric deformations,  as  in  simply-supported  and  clamped-clamped  beams,  because  the  resulting 
algebra  is  slightly  easier.  For  nonsymmetric  responses,  as  in  a sirnply-supported  clamped  beam, 
the  static  deformed  shape  should  be  used.  If  the  static  deformed  shape  is  not  used  in  non- 
symmetric cases,  uncertainty  will  otherwise  exist. 

In  all  of  the  illustrations,  the  assumed  deformed  shape  being  applied  to  a solution 
is  of  less  importance  to  the  resulting  deformations  and  strains  than  the  effects  of  coupling. 
Supporting  a llexible  structural  component  on  a flexible  foundation  has  a much  greater 
inllucnce  on  structural  response  than  the  assumed  deformed  shape,  as  the  next  part  of  this 
section  shows 

B.  Energy  Solutions  for  Coupled  Rigid-Plastic  Systems 


In  the  suppressive  structures  program,  energy  solutions  have  been  difficult  to  apply  be- 
cause the  actual  structure  is  a combination  of  plates  and  beams  or  a combination  of  I-beams 
and  hoops,  rather  than  a simple  beam,  plate,  or  membrane  configuration.  Although  energy 
solutions  developed  to  date  in  the  Suppressive  Structures  Program  (see  Refs.  1 , 4,  6,  and  9) 
apply  for  simple  structural  elements,  structural  configurations  which  are  combinations  of 
elements  can  also  be  solved  using  this  approach;  however,  to  complete  such  a solution,  one 
more  equation  is  needed  to  eouple  deformation  in  the  first  structural  element  with  defor- 
mation in  the  second.  For  a rigid-plastic  system,  derivation  of  the  required  relationship  is 
straightforward. 

ip  Coupled  plastic  systems  can  be  visualized  by  a rigid-plastic, 

rheological  model  as  in  Figure  .1.  In  Figure  1 , the  yield  force/i 
— I . . must  exceed  the  yield  force  /x , and  the  quasi-static  force  P must  be 
I m..  I — less  than  both/,  and /i : otherwise,  the  system  will  not  act  as  a 
I I coupled  system.  The  relative  residual  deformation  (A",  ~ Xi)  is 

I ^ related  to  the  residual  deformation  Xj  through  an  impedance  match 

f.  ||  involving  the  parameters , Wj./,  ./j,  and  F.  Equation  (34)  is 

* I 1 this  coupling  relationship. 


(-y.  --y.)  m,(f,-p) 

Xi  Wt  1 (/]  ■“  y 2 ) 

To  complete  an  energy  solution,  one  now  only  has  to  equate  the 
kinetic  energy  to  the  strain  energy  or  the  maximum  possible  work  to 
the  strain  energy,  to  obtain  estimates  for  the  deformation  in  a dynami- 
cally loaded  suppressive  shield.  In  a coupled  system  as  in  Figure  1,  the 
strain  energy  is  given  by: 

S.E.-/,(^,  -Afj )-!-/,  (.V,).  (35) 


i'lGURt  1.  RIGID-PLASTIC 
RIIUOLOCICAL  MODEL 


Because  two  equations  now  exist  for  and  Xj.  they  can  be  solved  simultaneously.  We  will 
proceed  by  developing  Eq.  (34)  and  illustrating  its  use  in  two  example  problems. 


19 


I.  Development  of  the  Coupling  Equation 


Assume  mt  is  hit  simultaneously  with  an  impulse  / and  a quasi-static  force  P.  Then 
the  two  equations  of  motion  are:* 


+/,=>'  (36a) 


WjAfj  4-/2  -ft 


(36b) 


For  no  initial  displacement,  no  initial  velocity  for  m2,  and  an  initial  velocity  of  I/m  i for  »i, , 
we  obtain  by  direct  integration: 


mt  2r/i,  2mi 
(fi 


(37a) 

(57b) 


Coupled  motion  continues  until  (x  1 — x 2)  equals  zero, 


(38) 


or  until  time  f ' given  by 


t' 


[(m,  +miy,  “m,/2-m2/’) 

At  this  time,  the  new  initial  conditions  for  m2  are: 

2[(m,  4- ffijyi -m,/2  — m2/’]’ 


(39) 


(40a) 


dxi 


(Ti  -h)l 


d/  [(w,  4- m2)/’i -m,/2  - W2^] 

The  relative  motion  (x,'  — x^)  equals  the  maximum  relative  residual  motion  (Xt  — X2) 

m2l^ 


(40b) 


(Xt  -X2)  = 


2m,  [(m,  4-  m2yi  “m,/,  ~ m2p] 


(41) 


After  element/,  “locks  up”  at  time  t' , motion  continues  as  an  uncoupled  system. 
The  equation  of  motion  in  this  second  phase  is  given  by: 


*Note  thut  / li  total  inipulw  with  dimeniloni  of  FT,  rather  than  specific  blast  impulic  f with  dimenilons  of  FT/L’. 


(m,  + = P 


(42) 


3y  establishing  u new  time  zero  ut  the  instant  l*  when  the  motion  becomes  uncoupled,  using 
the  initial  conditions  given  by  Eqs.  (40),  and  integrating,  wc  obtain  as  the  equation  for  dis- 
placement in  the  second  phase; 

^ -(y,  -P)t^  ^ u\  ~h)h  ^ t»,(/i  -A)/^ 

^ 2(«/|+Ht2)  1(»«1  + mj)/,  2[(/ni  -I- Wi)/!  -Wi/il* 

All  motion  stops  ut  time  rt  given  by; 


/t  = 


(W,  +W2)(/’|  -fi)l 

(/■j-/’)lvW,  +M7j)/i  — mi/j -mjP] 


(44) 


Substituting  /f  into  Eq.  (43)  and  gathering  terms  gives  the  maximum  residual  deformation, 
for  clement  mj. 


-f2V^ 

2{fi  -P)l(mi  + mi)fi  — ntifi  - ntiP] 


(45) 


Finally,  dividing  Eq.  (41 ) by  Eq.  (45)  gives  the  relative  displacement  coupling  equation 
already  presented  as  Eq.  (34). 

(X,  -Xi)  -P) 

X-i  I (/i  ~ ) 


This  equation  is  the  key  coupling  relationship  which  will  be  used  in  all  calculations.  Use  of 
the  equation  is  best  illustrated  by  sonnc  example  problems.  Calculation  of  the  mass  m for  a 
structural  component  is  no  problem;  however,  the  effective  force/ requires  assumptions  when 
calculated  for  a beam  or  plate  element.  We  will  calculate  the  force  /by  assuming  a deformed 
shape  for  a structural  component,  calculating  the  strain  energy  stored  in  that  deformed  shape, 
determining  the  average  deformation  for  a given  deformed  shape,  and  finally  dividing  the 
strain  energy  by  the  average  deformation.  Once  this  procedure  is  complete,  the  coupling 
Eq.  (34)  and  energy  relationships  permit  structural  deformations  and  strains  to  be  deter- 
mined in  a procedure  similar  to  that  used  in  any  uncoupled  structural  analysis. 


2.  Application  to  the  Category  I Shield 

The  Category  1 shield  is  used  as  the  approximate  illustrative  example  because  it 
was  the  first  to  draw  our  attention  to  the  need  for  a coupled  solution.  Here  we  will  derive 
the  equations  for  coupled  response  of  the  beams  and  rings.  In  Chapter  IV,  comparisons 
with  data  from  the  Category  1 l/4-sc3le  model  demonstrate  the  validity  of  the  coupled 
energy  solution. 

Basically,  this  shield  is  a barrel  as  shown  in  Figure  5 (page  52).  External  circular 
hoops  provide  restraint  for  longitudinal  I-beams.  The  I-beams  represent  the  m i — /,  structure, 
and  the  hoops  are  the  w,  - /j  structure.  We  will  treat  the  I-beams  as  claniped-clamped 
beams.  An  assumed  deformation  pattern  for  this  .structural  component  is  given  by; 
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off ' ' ;■  * 


I 


'I’ojiiiiitiffiiiwiipiiinfiiri 


(46) 


Wfl 


(l+c„s^) 


where 

Wo  = maximum  mid-span  deformation  relative  to  the  supports 
8 = total  span 

X ~ coordinate  system  with  origin  at  mid-span 
y = deformation  at  some  value  of  x 
Differentiating  Eq.  (46)  twice  gives: 

d^y  — 2ir^Wo  2ttx 

The  strain  Energy  S.E.q  is  given  by  the  integral 


8/2 


r d^y 

SE®  = 2 / My^^dx 


(47) 


(48) 


where 

My  - plastic  yield  moment 

Substituting  Eq.  (47)  into  Eq.  (48)  and  completing  the  desired  integration  give.s: 

The  average  deformation  ^avg  ^i^st  be  calculated  next.  It  is  obtained  by  the  integral 

8/2 


(49) 


I^IJ'avg 


KH  / 


2vrx  \ 

1 + cos  — 1 dx’ 


(50) 


or 


w„ 


Tavc  = 


(51) 


The  force /i  in  one  beam  equals  S.E.q  divided  by  the  average  deformation,  or: 


/.  = 


SirMy 


(52) 
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Tlic  mass  »i,  equals: 


/»!  = p*j4i,2 


(53) 


where 


Pb  ~ mass  density  of  a beam 
Ah  - cross-scctionai  area  of  a beam 


Next./  j and  «jj  must  be  calculated  for  the  hoops.  If  we  assume  hoops  of  equal 
effective  cross-sectional  areaM*,  on  each  end,  the  deformation  will  be  a symmetric  change 
in  radius  A/?.  The  strain  energy  for  both  hoops  is: 


S.E.q)  = (2ay ) 


(ItRAh) 


(54a) 


or 


S.E.^  = 47ra^/l;,  A/J  (54b) 

The  average  deflection  in  a hoop  is  the  deflection  A/?.  This  means  that  the  force /j  is  given 
by: 


A = S.E.@  /A/?  (55a) 

or 

A = ^TfOyAh  (55b) 

The  mass  Wj  equals  twice  the  mass  density  times  the  area  times  the  circumference,  or: 

wij  = (2pA)(>4/,)(27r)?)  (56a) 


or 


trii  - A-nphAhR  (56b) 

The  quasi-static  force  equals  the  pressure  p times  the  internal  circumference  times  the 
length,  or: 


/»=  (p)(27rr)(B)  (57a) 

p = 2jrrRp  (57b) 

Now  the  coupling  equation  can  be  used  to  relate  the  average  deformation  in  the 
beams  (rvu/2)  to  the  deformation  in  the  hoops  A/?.  Substituting  for  Wi , wij,  A >A  i snd 
P then  yields: 


2 i^TrphAhR)\{^TtahAh)-2vrip] 

(ZifNat>2  \ 

(ATpfcyiftB)  ^ — 4Tra,,Ahj 

The  parameter  N stands  for  the  number  of  beams.  The  mass  of  the  beams  and  force  equals 
iV  times  m,  or  J\ , respectively,  for  a single  beam  in  Eq.  (58).  In  addition,  the  yield  stress 
Ob  times  the  plastic  section  modulus  z was  substituted  for  the  yield  moment  in  the  beam 
force  equation.  Eq.  (52).  Reducing  Eq.  (58)  algebraically  gives  the  ratio  for  the  maximum 
deformation  in  the  beams  relative  to  the  change  in  hoop  radius. 


(59) 


The  solution  proceeds  by  writing  the  strain  energy  S.E.  for  the  entire  system. 
This  energy  equals  the  sum  of  Eqs.  (49)  and  (54),  or: 


S.E. 


AnNobZWo 

e 


+ AvohAbAR 


Substituting  Eq.  (59)  for  AR  in  Eq.  (60)  gives: 


(60) 


S.E. 


AnNobZWp  ^ NpbAbOh^Wo 

e 2phR 


(61) 


Next  the  energy  imparted  to  the  structure  must  be  estimated.  This  energy  comes 
from  the  kinetic  energy  imparted  through  blast  waves  and  the  work  from  the  quasi-static 
pressure  buildup  within  a suppressive  shield.  Algebraically,  this  energy  E\  equals  Eq.  (62) 
with  loads  imparted  to  the  beams. 


En  = + pt>i{yfcvc  + (62) 

where 

b = loaded  width  of  the  beams. 

Substituting  for  /,  m , , y avg  . and  Eq.  (59)  for  AR  yields,  after  collecting  terms: 


i 
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^ pbiWn  ^ ^ ^ NpbAtH  Vqft/<j,g  / 
2P6>»6  2 AnphAhR  I _P^\ 

\ lOhAh) 


(63) 


Finally,  equating  , Eq.  (63),  to  S.E.,  Eq.  (61),  and  gathering  terms  yields  an  equation 
whieh  can  be  solved  for  , the  maximum  deflection  in  the  beams. 


^ pb9}  ^ ^ ^ NPbAt^ 
PbAbOhZWu  Nat2  AvphAhR 


/2Nov  _ j \ 
Oh PbAb^^  Va^/l^g  / 

O TT  ' y r.  \ 

ObPh2R  /j__P£L) 

\ 2ar,Ah) 


(64) 

Equation  (64)  yields  w„,  the  maximum  beam  deformation  relative  to  the  rings, 
and  subsequent  substitution  in  Eq.  (59)  yields  A77,  the  change  in  hoop  radius.  Strains  can 
also  be  estimated.  The  residual  hoop  strain  equals  A/?//?  as  in  Eq.  (65a),  and  the  maximum 
residual  bending  strain  in  the  beams  equals  half  the  beam  depth  HI'l  times  the  maximum  beam 
curvature  as  given  by  Eq.  (47).  Equation  (65b)  is  the  maximum  residual  strain  in  the  beams. 


^hoop  ~ 

n^WoH 


^beuns 


(65a) 

(65b) 


Refer  to  Chapter  IV  for  calculations  based  on  these  equations  for  the  1 /4-scale  model  of  the 
Category  1 shield. 


3.  Application  to  the  Category  III  Shield 

The  second  example  will  be  that  of  a rectangular  membrane,  supported  rigidly 
along  two  opposite  edges  and  by  flexible  clamped-clamped  beams  along  the  other  two 
edges.  This  configuration  might  be  representative  of  the  original  Category  111  containment 
structure^'  which  was  replaced  by  the  1 /4-scale  Category  I shield.  We  begin  by  assuming 
a deformed  shape  for  the  membrane,  Eq.  (66). 


TTX  jr^ 

»v  = w„cos-cos- 


(66) 


In  a membrane  the  strain?  are  calculated  from  the  first  derivatives  of  the  slope  by: 


(67a) 


'■yy 


1 (hw\  ^ 1 f^x\  . ^/iry\ 

(ay/  " (2A'/®'"  [iy) 

(ow\  /dw\  Jr^vvj  /tix\  /jry\ 

srjW'TiT?”"!?)”"  [7) 


(67b) 


(67c) 


In  a structure  under  a biaxial  state  of  stress,  the  strain  energy  per  unit  volume  is  given  by 
the  integral: 


t68) 


If  we  assume  a rigid-plastic  material  with  a yield  point  in  the  plate  Op , then  Oxx  = Oyy  = Op 
and  Oxy  = Opl\/3  according  to  the  distortion  energy  yield  theoiy.  Substitution,  of  constant 
stresses  into  Eq.  (68),  integrating  Eq.  (68)  for  strains,  substituting  the  strains  from  Eqs.  (67) 
into  Eq.  (68),  and  expressing  the  volume  of  the  membrane  through  the  thickness  times  a 
double  integral  yields: 


/\  d\  \ / > / 


0 0 
2a ph  n^wl 

i6Ary 


,hTr^ 

SY^ 


(69) 


X cos* 


(5) -(5)] 


Performing  the  required  double  integration  gives: 


S.E., 


Tr^Ophwj  Y 


16  ;ir 

X ■^v^TT^  r 


"O  8 

The  average  plate  deformation  must  be  calculated  next  by  double  integration. 


(70) 


or: 


X Y 

//•  KX  Tty 

J dywo  cos  — cos—  (71) 

0 0 


4w„ 

»''avg  - f 


(72) 
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The  force /i  in  a plate  equals  S.E.0  divided  by  wavg  ■ 0‘  - 

_ T(*aphwo  fV  _16_  , ^ 1 

“ 32  LA'  Y J 


(73)- 


The  mass///,  equals: 


mt^ApphXY  (74) 

Next ,/j  and  ntj  must  be  calculated;  however,  this  has  already  been  done  in  the 
first  example.  Because  there  are  two  beams,  the  quantity equals  twice /in  Eq.  (52),  and 
the  mass  //ij  equals  twice  in  in  Eq.  (53).  These  results  are  suinmaiized  as  Eqs.  (75). 


l6ltMy 

s. 


t75a) 


/Ma  = 2pfc/4Q 


(75b) 


The  average  deflection  of  the  beams  is  still  m/a/2,  and  the  quasi-static  force  is  given  by  Eq.  (76). 


P=4pXY 


(76) 


We  arc  now  prepared  to  substitute  into  Eq.  (34).  Substituting  for  /n , , mj , /, , /j , and  P 
yields  the  coupling  equation  for  deflections. 


or; 


Su/p 


^IirphMyA  — SptipXYAS. 


^pphXY 


v*ar,hw„ 


32 


Y _16_  X / _ \6TtMy 
X ^ Y ^ 


(77) 


TT  ppOph^XYwl  \Y  'L  + - j - EpII^Y.^2. 
32  ptntyA  ]X^^^  ^ i 

pXYi  \ 

4My  ) 


(78) 


The  solution  proceeds  by  writing  the  strain  energy  S.E.  for  the  entire  system.  This  energy  is 
two  tiiries  Eq.  (49)  for  the  strain  energy  in  one  beam  plus  Eq.  (70). 


S.E. 


8irMy  Wb 
fi 


n^Ophwi  Y 16 
4- ? 1.  _ + 

8 X 


(79) 
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Substituting  Eq.  (78)  tor  und  collecting  terms  yields: 


.TT^gp/m-j  , 2p,_XYh 

8 PXYfl\  X Y 


[ISppMyXYhwo 

'--(-‘Sr) 

The  kinetic  energy  from  the  blast  wave  wlach  is  imparted  to  the  system  equals: 

iH^XY?  li^XY 
~ 2(4ppXYh)  pph 


The  work  from  the  quasi-static  pressure  loading  equals: 


WK=/)(4^1') 


Or  after  substituting  Eq.  (78)  for  w/, : 


\6pphXY 


Jr’p6>l£ 


(l-^ 

V 4Af* 


Equating  WK  + KE  to  S.E.  gives  the  final  solution.  In  nondimensional  format  this  solution 


IT/  + ITp  0.8106.7^,  + 


0.39277rp 


pVoVynl  _ 6484i^jv  _ q 
’To’Tp)  (4-»r<,jrp^ 


-\=  0.6169 


X 1 + 


(4  — iroir„) 


TTuTrJ,  --  81 .48 


riff  (4  irffTTp) 


PpXYh 


Oph^i 

PrXY 

“ Oph^ 

i^XY 

PpOph* 


16  X 
V'37r2  + Y 


(84b) 

(Cond) 


Equation  (84)  is  a quadratic  equation  that  must  be  solved  for  Jr.^  or  Wp . After  Wp  is  ob- 
tained, ivi,  is  obtained  from  Eq.  (78),  which  in  nondimensional  format  is: 


Wp  _ Q.^92lKg‘n„tyitv,  — 6.484rrp 


(4-ir(,irp) 


(85) 


To  comolete  this  illustration,  let  us  apply  Eqs.  (84)  and  (85)  to  the  Category  III 
shield.  We  will  assume  parameters  as  follows; 

p = 66  psi 

i * 0.179  psi-sec 

X = 26.5  in. 

y=  56.5  in. 

h = 1 .033  in.  (this  is  an  effective  thickness) 

. lb-sec^ 

Pp--=Pb  = 7.33  X 10-“  — T- 


Ob  - 48,000  psi 
A = 61.71  in.’ 
fi  = 113  in. 

My  = 1.21  X 10*’  in.-lb 
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Substitution  of  these  parameters  into  the  it  terms  gives: 


Up 

= 0.2218 

T^o 

= 0.4783 

Itp 

= 1.929 

Vi 

= 1.197 

Ity 

= 3.537 

Substituting  into  Eq.  (84)  gives  a positive  root  for  riw  of  3.953;  hence,  Wp  = 4.083  inches. 
Finally,  substitution  into  Eq.  (85)  gives  wt  - ~1  135  inches.  This  negative  sign  associated 
vvith  wi,  is  not  a mistake.  It  means  that  the  plate  arid  beam  systems  are  not  coupled  as  has 
been  assumed;  hence,  this  analysis  is  not  appropriate  for  the  parameters  substituted  into  it. 
In  this  particula,  solution.  /]  is  greater  than/i . This  means  there  is  no  plastic  deformation 
in  the  second  element,  and  the  particular  system  under  investigation  is  uncoupled.  Appro- 
priate answers  can  be  obtained  by  using  uncoupled  techniques  which  have  already  been 
presented.  The  purpose  of  this  second  illustration  was  to  emphasize  the  meaning  of  a 
negative  sign  when  it  is  encountered  in  this  type  of  analysis. 

4.  Importance  of  Treating  Coupled  Response 

We  have  shown  that  enelgy  solutions  are  possible  for  coupled  perfectly  plastic 
systems.  For  two  members  in  a coiSpled  system,  one  additional  equation  is  needed  to  relate 
deflections  in  the  first  member  to  deflections  in  the  second  member.  Equation  (34)  is 
presented  as  a relationship  for  coupling  the  deflections.  With  the  assistance  of  this  one 
extra  equation,  energy  solutions  can  be  developed  as  in  any  rigid-plastic  coupled  system. 

I 

In  Chapter  IV  wc  will  show  that  the  effect  of  coupling  is  very  important  for  pre- 
dicting response  in  the  Category  I shield.  On  the  other  hand,  the  example  above  shows 
that  coupling  is  unimportant  for  predicting  membrane  deformations  in  the  Category  III 
shield.  While  these  solutions  are  more  complicated  than  uncoupled  solutions,  residual 
deilec  Lions  and  strains  for  coupled  systems  can  be  predicted  with  closed-form  algebraic 
solutions,  and  the  usefulness  of  the  approximate  energy  methods  for  predicting  response 
to  blast  or  impact  loading  is  greatly  enhanced. 

C.  Elastic-Plastic  Energy  Solutions  for  Beams 

In  the  design  of  suppressive  shields  we  are  usually  interested  in  predicting  onset  of 
failure,  and  hence  maximum  strain  or  maximum  stress.  In  most  problems,  a beam  has 
segments  which  remain  elastic  and  segments  which  go  plastic.  The  assumptions  of  total 
elasticity  or  total  plasticity  throughout  all  elements  in  a structural  component,  or  concen- 
trated plasticity  in  a hinge,  are  approximations.  To  obtain  an  idea  of  how  elastic-plastic 
solutions  could  be  presented  to  designers,  and  to  see  how  combined  solutions  might  be 
obtained  from  simple  elastic  and/or  simple  plastic  solutions,  two  elastic-plastic  solutions 
have  been  developed  for  a simply-supported  beam  and  a string. 
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The  constitutive  relationsltlp  used  Cor  both  solutions  is  given  by  Eq.  (86). 

o-  Ov  tanh  — . (86 

■ 

This  is  an  excellent  approximation  to  an  elastic-plastic  malerial.  For  small  strains,  Ee/oy  < 
0.5.  the  stress-strain  curve  is  a linear  one  with  o = Ee:  and,  for  large  strains,  Eeloy  > 2.0, 
the  stress-strain  curve  is  a perfectly  plastic  one.  with  a Oy. 


1.  Solution  for  a Simply-Supported  Beam 

in  any  structural  element  under  a uniaxial  state  of  stress,  the  strain  energy  per 
unit  volume  is  given  by: 


= foi 


If  we  substitute  Eq.  (86)  for  a and  integrate  from  a strain  of  zero  to  the  maximum  strain  e, 
we  obtain; 


--- 


U _Oy^ 
Vol.~  E 


in  cosh I 


Equation  (89)  is  then  applied  to  any  bending  or  extensional  solution  by  substituting  the 
strain  and  integrating  over  the  volume.  If  we  use  the  assumed  deformed  shape  for  a simply- 
supported  beam  given  by  the  first  term  of  Eq.  (1 ),  the  strain  e for  a differential  element  in 
bending  is: 

, d*y  Tt'^Woh  Ttx 

€ = —h  — sin  (90 


If  the  beam  is  of  rectangular  cross  section,  the  volume  equals: 


Vol.  = 4b 


Un  8/2 

/ d*  / 


Substituting  the  expression  for  e and  the  volume  integral  into  Eq.  (89)  gives  the  following 
double  integration  for  the  strain  energy  £/.  . . 


4ai,^b 


Hl2  8/2 

/ f 


in  cosh 


(v^WoEh  . 7rx\ 
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or 


<7  = 


4ajfe//g 

7r^£ 


(93a) 


0 0 


where 


a = TTjf/g  and  0 ® ith/H. 


(93b) 


Next  we  assume  that  the  loading  is  a uniformly  applied  impulse,  giving  for  the  kinetic 
energy: 


2pH 


(94) 


Equating  the  kinetic  energy  to  the  strain  energy  and  rearranging  yields  the  dimensionless 
equation: 


where 


(95) 


(96a) 


C~{it^(^El%paylP)9.nAD  = {-irw„HEIay9})  (96b) 

A computer  program  was  required  to  numerically  perform  the  double  integration.  Results 
of  this  program  are  represented  by  the  solid  continuous  line  in  Figure  2.  Dashed  jines  in 
the  figure  are  the  asymptotes  for  completely  elastic  or  fully  plastic  behavior.  Note  that  the 
elastic-plastic  solution  correctly  approaches  the  elastic  asymptote  for  small  deformations 
and  the  plastic  asymptote  fo.'-  large  deformations. 

2.  String  Solution 

The  string  solution  uses  the  same  strain  energy  per  unit  volume  equation,  Eq.  (89), 
as  the  beam  solution.  These  solutions  differ  in  the  equations  relating  strain  to  deformation. 
In  an  extensional  element  the  strain  is  given  by  Eq.  (20a).  Substituting  the  first  term  of 
Eq.  ( 1)  for  the  assumed  deformed  shape  in  Eq.  (20a)  gives  Eq.  (97). 
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Tlie  volume  for  un  extensionul  element  equalsi 

Vol.  = 2A 

where  A is  the  cross-seetiunal  area. 


e/3 

y dv 


(98) 


Substitutine  Kqs.  (97)  iind  (98)  into  Eq.  (89)  gives  the  following  integral  equation  for  the 
strain  energy  U. 


jr  to 


(99) 


or.  using  Eq.  (93b), 


“1 

j in  cosh  ^ 282  0^ 

j cos*  aj 

(100) 


Thckinetic  energy  imparted  to  a string  equals; 


2rtt  2pA 


(101) 


Equating  the  kinetic  energy  to  the  strain  energy  and  rearranging  yields  the  nondimcnsional 
equation; 


(102) 


This  equation  can  be  written  as: 


F=  I fin  cosh  [ O cos*  a]  da 

•o 


(103a) 


where 


F=  (nF t'/Apa^A^ ) and  G = (jr*w3£'/2S*ay)  (103b) 

Numerical  integration  of  Eq.  (103)  yields  the  solid  line  in  Figure  3.  As  for  the  beam  solution, 
the  dashed  lines  in  Figure  3 represent  the  asymptotes  for  fully  elastic  or  fully  plastic  behavior, 
and  the  inset  is  a nondimcnsional  plot  of  Eq.  (86). 
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FIGURE  3.  ELASTIC^LASTIC  SOLUTION  FOR  A STRING  LOADED  BY  A UNIFORM  IMPULSE 


Although  the  elastic-plastic  solutions  given  in  Figures  2 and  3 are  extremely  useful 
for  designers,  they  are  somewhat  difficult  to  derive  because  numerical  integration  is  required. 
We  will  now  show  that  the  numerical  solution  can  be  approximated  to  within  one  or  two  per- 
cent by  properly  combining  the  elastic  and  plastic  asymptotes.  The  asymptotes  are  easily 
derived. 


.1  Limiting  Elastic  and  Plastic  Cases 


Before  an  approximate  elastic-plastic  solution  can  be  found,  the  asymptotes  are 
required.  Although  these  limiting  cases  have  been  derived  many  times  before  (see  Refs.  1 . 
4,  6,  and  9),  we  will  rederive  them  here  for  completeness. 


In  elastic  problems,  the  strain  energy  per  unit  volume,  Eq.  (87)  becomes: 


JL-i.  i 

Vol."  2 ^ 


(104) 


In  plastic  problems,  the  strain  energy  per  unit  volume,  Eq.  (87)  becomes: 


U 

Vol. 


= OyC 


(105) 


Of  course,  these  equations  apply  to  either  the  beam  or  string  solution. 


For  the  beam  solution,  both  elastic  and  plastic  strains  are  given  by  Eq.  (90), 
and  the  volume  equals  the  double  integral  given  by  Eq.  (91 ).  Substituting  these  two  equations 
into  Eq.  (104)  for  the  elastic  case  or  Eq.  (lOS)  for  the  plastic  case  give«; 


elastic: 


U = 


2r,*Ew„-^b  "c  V 


/ / sin» 

0 0 ' / 


dh  dx 


rigid-plastic; 


tf/2  B/2 


„ 4ir»a,w„b  T ; . . /^A  ..  . 

^ 8^ J J /rsinfjjdAtU 

r\  ' ' 


0 0 


Or,  after  completing  the  required  double  integration:;: 


(106) 


(107) 


elastic: 


C/  = 


48fi* 


(108) 
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rigid-plastic: 


noy  Wq bH^ 


6'  = 


2K 


(109) 


Now,  by  equating  Eq.  (108)  to  Eq.  (94)  and  Eq.  ( 109)  to  Eq.  (94).  the  elastic  and  plastic 
asymptotes,  respectively,  of  Figure  2 are  obtained.  These  are  given  by  Eqs.  (110)  and  (111). 


(elastic  beam) 


tH*' 

pEwo^H*  24 


(110) 


and 


(rigid-plastic  beam) 


i‘^9} 

POyWoH^ 


(111) 


For  the  string  solution,  both  elastic  and  plastic  strains  are  given  by  Eq.  (97),  and 
the  volume  is  given  by  the  integral  in  Eq.  (98).  Substituting  these  two  equations  into  Eq. 
(104)  for  the  elastic  string  or  Eq.  (105)  for  the  plastic  string  gives; 


elastic: 


48“ 


(112) 


rigid-plastic: 


ir^OyWo^A 

8^ 


/ 


Or,  after  completing  the  required  integrations; 


elastic: 


(/  = 


3k*EAwo* 

648’ 


(113) 


(114) 


37 


rigid-plustic: 


£/=« 


n^OyAwp^ 

4C 


(115) 


The  elastic  and  plastic  asymptotes  for  the  string  solution  are  obtained  by  equating  Eqs.  (114) 
and  (115),  respectively,  to  the  kinetic  energy  given  by  Eq.  (101).  Equations  (1 16)  and  (117) 
obtained  in  this  way  are  the  asymptotes  of  Figure  3. 


(elastic  string) 


fibril*  3»r* 
pEA^Wp*  ~ 32 


(116) 


and 


(rigid-plastic  string) 


pa^>4*Wo*  2 

The  approximate  elastic-plastic  solutions  can  now  be  derived. 


(117) 


4 Approximate  Elastic-Plastic  Solutions 

The  elastic-plastic  solutions  in  Figures  2 and  3 are  very  closely  approximated  by 
Eq.  ( 1 1 8)  where  Y is  the  ordinate,  X is  the  abscissa,  and  coefficients  A and  B are  constants. 

Y ^ A Unb}  (BXy  ' ^ (118) 


The  elastic  asymptote  for  Eq.  (118)  is: 

(elastic  asymptote) 

J = (^5)  (119a) 

And  the  plastic  asymptote  is : 

(plastic  asymptote) 


y=(/4)  (119b) 

These  equations  have  the  same  form  as  the  asymptotes  for  the  beam  and  the  string.  Thus, 
for  the  simply  supported  beam  the  terms  in  Eqs.  (119)  are  defined  by  comparison  with 
Eqs.  (1 10)  and  (1 1 1).  We  find: 
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y s 

POyWoH^ 


(120a) 


IH* 

" pEwo 

(120h) 

/4  = ff 

(120c) 

ll20d) 

Solving  for  X and  B and  substituting  into  Eq.  (118)  gives  an  accurate  approximation  to  the 
elastic-plastic  bending  solution  without  encountering  the  inconvenience  of  a numerical 
integration  of  a complicated  function. 


pOyWoH^ 


= 3.1416  tanh* 


(121) 


To  compare  with  Figure  2,  this  solution  can  be  recast  in  the  format: 


vEwoH\  , 


(122a) 


or 

C=  1.234Z)tanh*  (0.64127*'*)  (122b) 


Either  of  the  preceeding  formats  is  acceptable  for  the  elastic-plastic  beam  solution. 
The  advantage  to  the  format  given  by  Eq.  ( 1 22)  is  that  the  deformation  Wg  is  isolated  on  one 
side  rather  than  appearing  on  both  sides.  Designers  would  prefer  to  isolate  Wg  in  this  manner 
because  this  deformation  also  directly  determines  the  maximum  strain.  Equation  (122)  is  not 
the  computed  line  appearing  in  Figure  2;  however,  it  is  very  difficult  to  detect  differences  in 
Eq.  (122)  and  the  computer  solution. 

Equation  (118)  and  its  asymptotes  can  also  be  used  to  obtain  an  accurate  approxi- 
mation to  the  string  solution.  The  elastic  and  plastic  string  solutions  are  given  by  Eqs.  (116) 
and  (117).  Comparing  Eqs.  (1 19)  to  these  equations,  we  have: 


Y (*fr*g* 
x"  pEA^Wo^ 

(*6*8* 


y = 


POyA^Wo^ 
3ir* 


AB  = 


32 


(123a) 

(123b) 

(123c) 
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{123d) 


Solving  forfl  and  X and  substituting  into  Eq.  (118)  yields; 


i^bH^  ir^  

-j — : = — tanh*  2.356  ( — ) 

pOyA^Wo^  2 yay9}  j 


(124) 


Or  after  modifying  the  format  to  isolate  terms  containing  rvo  on  one  side  and  using  Eq.  (103b), 
'■ni^Eb^ 


Apay‘^A'‘ 


= 0.7854 


ir^Ewo^~\  f fir^Ewn^\ 

^Ju„h>  ^0.6,34  J ““•) 


or 


F=  0.7854  G tanh^  (0  6 124  G''^) 


(125b) 


This  solution  also  approximates  the  more  detailed  elastic-plastic  solution  with  the 
same  degree  of  accuracy  as  the  bending  beam  analysis.  Philosophically,  one  should  not  be 
disturbed  by  the  use  of  the  hyperbolic  tangent  squared  as  a function  for  combining  elastic 
and  plastic  solutions.  After  all,  any  solution  requires  an  approximate  stress-strain  curve. 

No  approximation  to  a stress-strain  curve  uniformly  matches  all  materials.  Instead  of 
approximating  the  stress-strain  curve,  one  can  elect  to  approximate  the  strain  energy  trans- 
ition from  elastic  to  plastic  as  we  have  done. 


D.  Graphical  Solutions  for  Beams 

The  procedures  developed  in  the  preceding  section  can  be  used  to  derive  general  graphical 
solutions  for  blast  loaded  structural  elements.  These  graphs  are  attractive  because  they  permit 
rapid  solutions  to  difficult  problems  without  recourse  to  complicated  mathematical  procedures 
In  fact,  as  we  will  demonstrate,  the  graphical  solutions  are  self-contained  and  can  be  easily 
applied. 

Our  example  for  these  solutions  will  be  beams  loaded  in  the  same  manner  as  those  in 
the  case  of  the  Category  1 shield.  These  beams  are  subjected  to  an  initial  impulse  produced 
by  the  shock  wave,  plus  a long  duration  quasi-static  pressure  produced  principaily  by  the 
heating  of  air  in  a confined  space.  Because  the  blast  wave  part  of  the  loading  is  described 
only  by  an  initial  impulse  (not  a pressure  and  impulse  or  pressure  and  time),  this  solution  is 
suitable  only  when  the  duration  of  the  overpressure  in  the  blast  wave  is  less  than  one-quarter 
of  the  fundamental  period  of  the  beams.  The  pressure,  p,  in  this  solution  refers  to  the  quasi- 
static  pressure,  not  the  blast  wave  overpressure. 

In  response  to  these  loads,  the  beams  undergo  elastic  and  plastic  deformations.  For 
this  example,  coupling  between  the  beams  and  rings  is  neglected,  but  different  end  conditions, 
i.e.,  simply-supported,  clamped  or  a combination  of  the  two,  are  included. 
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The  soluiion  to  this  problem  is  similar  to  the  impulsively-loaded  beam  solved  in  Section 
C except  that  an  additional  term  is  required  to  define  the  energy  produced  by  the  action  of 
the  quasi-static  pressure,  p,  on  the  beam  as  it  displaces.  With  the  addition  of  this  term  called 
the  work.  WK,  the  energy  balance  is  written  as 


Kh + h‘K--U  (126) 

KE  is  the  kinetic  energy  produced  by  the  initial  impulse,  and  U is  the  strain  energy  (elastic 
and  plastic)  stored  or  absorbed  by  the  beam. 

Kinetic  energy  is  the  same  for  each  beam  regardless  of  the  support  condition  and  is 
given  by  Eq.  (94),  derived  previously.  The  work  and  kinetic  energy  depend  upon  the  mode 
shapes  assumed,  and  this  will  differ  for  each  support  condition.  For  these  calculations,  we 
have  assumed  that  the  deformed  shape  under  dynamic  loading  is  the  same  as  the  static 
deformed  shape.  The  calculations  are  illustrated  in  detail  for  a simply-supported  beam. 

The  static  deformed  shape  for  a uniformly-loaded,  simply-supported  beam  is  given  by 
Eq.  (127) 


y = ^^(x‘' -2Cx3  + H^x)  (127) 

where  iv„  is  the  center  deflection  (a-  = E/2).  We  tan  now  derive  the  work  from  the  quasi- 
static pressure  p as 


WK  = 


8 

J*  pby  dA 


(128) 


where  b is  the  beam  loaded  width  and  E is  its  length.  Substituting  Eq.  (127)  into  (128)  and 
performing  the  required  integration,  the  work  is  found  to  be 

= (129) 

As  shown  in  Eq.  (90),  the  strain  in  a beam  cross-section  at  a distance  h above  or  below 
the  neutral  axis  is 


€-—h 


d^y 

dx’ 


Again  using  the  deformed  shape  of  Eq.  (127),  we  have 


€ 


-h 


16wo 

5E* 


(i2A*  - 12Ex) 


(130) 


(131) 
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Thu  strain  unurgy  cun  now  be  computed  from  Eq.  (89);  however,  the  solution  can  be  made  more 
genera!  if  we  rel'orinulate  the  strain  energs'  in  terms  of  gross  cross  section  properties  of  the 
beam.  This  can  be  don**  by  replacing  the  stress-strain  constitutive  relationship  given  in  Eq.  (86) 
by  an  approximation  fur  the  moment-curvature  relationship.  Such  a relationship  is  given  by 
Eq. (132) 

£/v" 

yW  ^ A/p  tanh  — 7-  (132) 

Mp 

where  j " is  the  second  derivative  of  the  beam  displacement  with  respect  to  x (curvature),  and 
M,,  is  the  fully  plastic  moment  of  the  beams.  Notice  that  Eq.  (132)  has  exactly  the  form  of 
Eq.  (82).  It  was  not  derived  from  Eq.  (82)  and  represents  a slightly  different  appro.ximation 
to  the  stress-strain  behavior.  In  fact,  the  stress-strain  behavior  will  differ  slightly  for  beams 
of  different  cross-sections. 

The  strain  energy  per  unit  length  for  a beam  in  bending  can  be  written  as 


u c /£iy"\  /Ely"  \ 

- j Mp  tanh  { — — ) d.v  = — f fin  cosh  ( — - ) 

dx  ” V } E!  \Mp  ) 


(133) 


For  the  total  strain  energy,  Eq.  (133)  i.s  integrated  over  the  beam  length,  fi,  to  obtain: 


U-  f/l  tocosh  d.. 


Substituting  for  .v"  from  Eq.  (127)  gives 


(134) 


El 


V 

r 

(fJwX  192  , 

1 cosh 

0 

dx 


(135) 


It  is  convenient  to  nondimensionalize  the  integral  by  letting 


^ dx-  = Kdt 


(136) 


With  these  substitutions  wc  have 


(/  = 


A/gfi 

£■/ 


(’.37) 
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Now  substituting  Eqs.  (94),  (129),  and  (137)  into  Eq.  (126),  we  obtain  Eq.  (138).  which 
relates  the  deformation  of  a simply^supported  bean!  to  its  basic  properties  and  the  blast 
loading  parameters: 


16  _ Mji  r. 

+ f'biiwo  = -ff  J ««cosh 


2pA  25 ' 


KEIwo\  192 
MpeV  5 


te  -5) 


dt 


(138) 


Eo.uation  (138)  is  nondimensionalized  by  dividing  each  side  by  the  coefficient  of  the  integral, 
or: 


^ \6pbElWo 


2M}A  ISM} 

In  nondimensional  pi  terms,  the  equation  is: 


u,  + = y in  cosh  (t*  ‘"tjj  <^5 


(140) 


A more  convenient  relationship  is  obtained  if  w©  is  eliminated  from  all  but  one  group.  Tnis 
is  accomplished  by  rewriting  ttj  as  the  product  of  and  a new  group  ir'j . 


, {\6pbfl^\  fEJwA 


(14i) 


Also,  if  one  wishes  to  limit  the  maximum  strain  in  the  beam  rather  than  the  maximum 
deformation,  Wq  can  be  replaced  by  fmax  in  using  Eq.  (j31).  For  a simply-supported 
beam,  occurs  at  x = ill  and  h - Hfl.  Thus, 


With  this  substituticn,  becomes 


Now  Eq.  (139)  becomes: 


TTj  = 


24Hwo 


5£/e„ 


(142) 


24MpH 


(143) 


1 /Ibj/ET  y ^ /Pbi^  \ (EIe^\  r 

2\Mpy/^/  \apMp/  \<x^Mph)  J 


(144) 
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where  the  constants  in  rri  and  ^3  have  been  replaced  by  Op  and  a,,  respectively. 

Similar  equations  can  also  be  derived  for  clamped-clamped  and  clamped-simply-supported 
beams.  If  wc  again  use  the  static  deformed  shape  under  uniform  loading,  the  following 
equations  are  obtained. 

Deformed  Shape 


c-c: 

c-ss:  y - -4  £x»  + i £^x) 

ic*’  ^2 


Strain  Equation 

32w 

c-c:  e = — ° (6x*  — 6fijc  + 8^ ) 
7 7w 

c-ss:  e = — — 9£)f) 


Maximum  Strain 


c-c:  Cm*)! 


9} 


ll.SSHwo 
c-ss:  Cmax  “ g2 

The  defoi'med  shapes  are  substituted  into  Eq.  (128)  to  compute  the  work,  and  the  strain 
expressions  are  substituted  into  Eq.  (134)  to  compute  the  strain  energy.  These  expressions 
are  then  combineu  with  the  kinetic-energy  (Eq.  (94)]  according  to  Eq.  (126)  to  obtain  ex- 
pressions similar  to  Eq.  (138).  Performing  the  manipulations  described  for  simply-supported 
beams  and  generalizing  the  results,  on  equation  can  be  obtained  which  applies  to  all  three 
boundary  conditions.  It  is  given  by  Eq.  (145). 

1 / lb^\  ^ ^ jpbi^  \ I EJe^,^  \ _ 

2 \Mp^J  \apMp)  \a,MpH  ) ~ 

X(C,{*  +C,t  + C,)J  dt 

Note  that  the  only  difference  between  Eqs.  (145)  and  (144)  is  in  the  description  of  the 
nondimensional  deformed  shape.  The  proper  equation  for  each  beam  boundary  condition 
is  obtained  by  substitution  of  the  appropriate  constants  aii  defined  in  Table  111; 


f in  cosh  1^ 


EIe„ 


a,MpH 
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TABLE  III.  CONSTANTS  FOR  EQUATION  c 141) 


“n 

o* 

c^ 

c, 

1 

damped  (e-c) 

1.875 

16 

192 

-192 

+32  ‘ 

dampcd-simply  supported  (c-ss): 

1.732 

11.55 

92.4 

- 63.3 

0 i 

simply  supported  (ss-ss): 

1.5625 

4,8 

38.4 

- 38.4 

0 

i 


Equation  (145J  was  numerically  integrated  for  each  beam  boundary  condition  to  obtain 
the  graphical  solution  of  Figure  4.  To  better  understand  the  information  presented  in  the 
figure,  rewrite  Eq.  (145)  in  the  following  form; 

rr^  + tt'iTTj  = I (146) 

where  I is  the  value  of  the  integral.  Note  that  if  the  applied  pressure,  p,  is  zero,  then  the 
dimensionless  group  tt'j  is  zero,  and  Eq.  (146)  becomes 


]_ 

2 


rr?  =/ 


(147) 


This  equation  yields  the  asymptote  for  impulsive  loading,  or 


^i=ff||pao=v^  (148) 

Likewise,  for  ( = 0,  tt,  =0,  and  the  asymptote  for  pressure  loading  becomes 


(149) 


Equation  (146)  can  be  written  in  terms  of  these  asymptotes  if  it  is  multiplied  by  1//.  This 
gives 


(150a) 


(150b) 


(150c) 


This  is  the  general  form  plotted  in  Figure  4.  Because  the  integral  / depends  upon  tts  and 
the  beam  boundary  conditions,  the  asymptotes  .4 , and  /Ij  given  in  the  inset  as  a function 
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FIGURE  4.  D-SIGN  CHART  FOR  ELASTIC-PLASTIC  BEAMS  SUBJECTED  TO 
INITIAL  IMPULSE  PLUS  A QUASI-STATIC  PRESSURE 
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of  TT]  und  the  beum  boundary  conditions.  One  difference  in  the  graphical  solution  and  the 
form  given  by  Eq.  (145)  is  the  absence  of  the  constants  Up  and  a,.  Because  separate  curves 
are  plotted  for  the  different  boundary  conditions,  these  constants  could  be  eliminated.  For 
example,  to  eliminate  Op,  the  rri  asymptote,  ,4, , was  simply  multiplied  by  the  constant 
given  for  Up  in  Table  ill  for  the  appropriate  boundary  condition. 

Observe  that  the  solution  for  either  impulsive  loading  only  or  for  pressure  loading  only  is 
given  directly.  This  is  clear  from  Eq.  (150b)  and  from  Figure  4.  If  p “ 0,  for  example,  then 
TTj  = 0 and  we  have 


^ (151 

Thus,  if  the  beam  properties  and  the  impulse  are  known,  can  be  computed  and  xfj  read 
directly  from  the  appropriate  curve  in  Figure  4.  Alternately,  if  one  wished  to  know  what 
impulse  would  produce  a prescribed  strain  in  a beam,  ttj  could  be  evaluated,  i4 , read  from 
the  graph,  and  I computed  from  Eq.  (151). 


A similar  approach  holds  for  pressure  loading  only;  however,  for  the  case  of  an  initial 
impulse  plus  a quasi-static  pressure,  an  iterative  approach  is  required.  Consider  as  an  example 
the  beams  in  the  cage  of  the  Category  I 1 /4-scale  model.  Beam  properties  and  the  loading 
parameters  are  listed  below.  (Refer  also  to  Chapter  IV.) 


Beam  Properties: 

A~  1.67  in.^ 

/=  2.52in.'» 

Zp  - 1.95  in.'* 
b~  1,43  in.  (loaded  width) 
//•=3.0iii. 

£-30  in. 

Oy  = 45,000  psi 
30  X 10®  psi 

p = 0.000733  lb-sec’ /in. ^ 


With  these  values,  tt,  and  rr'i  can  be  computed. 


Loading: 


i " 0.48  psi-sec 
p = 190  psi 


ibyJWr 

TT,  1,944 

f’yZps/pA 

, pbi^ 

--=2.787 

OyZp 

To  find  the  maximum  strain  in  the  beams,  first  compute  the  values  for  pressure  and  impulse 
separately. 

Impulse  only:  /t)  =ir|  = 1.944 
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From  Figure  4 for  n clamped  (c-c)  beam: 


Jfj  = — 7;*3.2 

Oy  ZpM 


^max  ” 0.01  114 


Pressure  only;  /4j  - tt'j  = 2.787 


From  Figure  4, 


0.235 

OyZpH 


fimax  = 0.000818 


The  beam  is  more  sensitive  to  the  impulse  than  the  pressure.  Also,  the  maximum  strain  for 
the  combined  loading  will  be  larger  than  the  sum  of  the  strains  computed  separately.  There- 
fore, try 


For  this  strain,  itj  is; 


and  from  Figure  4, 


fmax  =1.1(0.01114  + 0.0008 1 8)  = 0.0 1 1 96 


rr,  = 3.434 


A,  =2.02  /li=17.9 


A I and  /4i  must  satisfy  Eq.  (150b)  (or  the  curve  in  Figure  4),  which  for  this  problem  becomes 


Substituting  A , and  A 2 gives 


For  a closer  result,  try 


/l.944\  ^ /2.787  \ 

iTfrj  *v7r)‘' 

i 

/l.944\=‘  /2.787  \ 


1.944 

V*‘°  17.9 


I 


From  Figure  4,  the  corresponding  values  of  >4  j and  ttj  are  . 

<42=18.1 
JTj  = 3.65 

Now  check  Eq.  ( 1 50b)  again. 

/1.944\  » 2.787  . . 

Urj 

This  result  is  close  enough.  Therefore,  the  maximum  strain  for  the  S3  X 5.7  beams  in  the 
Category  I 1/4-scule  model  for  clamped  ends  and  rigid  support  is  predicted  to  be 

emax  = 3.65  - 0.01271  = 12,710^6 

Also  note  that,  with  the  asymptotes /I  | and  deteimined,  the  entire  p-i  diagram,  for  a 
maximum  strain  of  12,710/re  in  the  beam,  can  be  established  from  Eq.  (150b)o;  from  Fig- 
ure 4.  So,  what  we  have  really  obtained  is  a general  solution  for  the  beam  and  not  just  the 
solution  for  one  set  of  loading  parameters. 


If  one  wishes  to  compute  maximum  deflection  in  the  beams,  the  relationship  between 
maximum  strain  and  the  deflection  can  be  used.  For  a simply-supported  beam,  this  relation- 
ship is  given  by  Eq.  (142),  and,  for  the  other  boundary  conditions,  by  the  expressions  on 
page  44.  Generalizing,  we  can  write 


Wo 


a,// 


(1S2) 


where  a,  is  given  in  Table  III. 

Thus,  we  have  obtained  a general  graphical  solution  which  includes: 

• beams  of  different  boundary  conditions  (with  rigid  support) 

• elastic-plastic  beam  behavior 

• predictions  of  maximum  strains  and  deflections 

• beam  response  predictions  for  any  combination  of  an  initial  impulse,  /,  and 
quasi-static  pressure,/?. 

Such  a solution  provides  a convenient  tool  for  preliminary  design  of  beams  subjected  to  this 
particular  type  of  loading.  For  final  design,  a more  rigorous  approach  is  usually  required, 
which  inclu  des  the  effects  of  support  tlexibility,  etc.  The  need  for  a more  rigorous  analysis 
will  be  demonstrated  by  the  comparisons  in  Chapter  IV  between  measured  and  predicted 
strains  for  the  Category  1 1,4-scale  model. 
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III.  EQUATIONS  FOR  RESPONSE  OF  STRUCTURAL 
ELEMENTS  TO  BLAST  LOADING 


Numerous  equations  for  estimating  the  response  of  structural  components  to  blast 
loading  have  been  developed  during  the  Suppressive  Structures  Program.  Derivations  of 
these  equations  have  been  reported  by  Baker,  et  al..^'  ’ by  Westine  and  Baker,^**  and  by 
Westine  and  Cox.'"*’  A partial  summary  was  included  in  the  report  by  Baker,  et  al.^®* 

Here,  a final  summary  is  made  which  collects  the  equations  reported  in  References  1, 4, 
and  6,  plus  other  equations  that  have  not  previously  been  reported.  Equations  presented 
are  for  the  response  of  individual  structural  elements.  Rigid  behavior  of  the  component 
which  loads  the  element  or  the  support  for  the  clement  is  assumed.  Equations  for  coupled 
responae  have  been  developed  for  a few  cases,  but  these  are  included  separately  in  Chapter 
II.  Chapter  II  also  includes  graphical  solutions  for  selected  components 

This  summary,  given  in  Table  IV,  is  an  expansion  of  Table  B-1  in  Reference  9.  Cor- 
rections to  equations  have  been  made,  as  required,  and  both  elastic,  elastic-plastic  and  rigid- 
plastic  solutions  are  included.  Structural  elements  covered  include  beams,  rings,  membranes, 
plates,  cylinders,  and  spheres. 

Equations  in  Table  IV  relate  the  peak  deformation  in  the  structural  element,  usually 
designated  as  , to  the  element’s  material  and  geometric  properties  and  to  the  applied  loading. 
The  applied  loading  is  treated  as  the  simultaneous  application  of  an  impulse,  ir.  and  a quasi- 
static pressure,  p.  This  loading  is  representative  of  that  on  suppressive  structures  produced 
by  the  internal  detonation  of  a high  explosive.  The  blast  wave  from  the  detonation  is  of 
short  duration  relative  to  structural  frequencies  (and  so  can  be  treated  as  an  impulse),  and 
heating  of  the  air  in  the  structure  produces  a pressure  buildup  of  much  longer  duration. 
Simultaneous  application  of  the  impulse  and  quasi-static  pressure  is  supported  by  pressure 
measurement  from  the  I / 1 6- scale  venting  tests  reported  by  Schumacher  and  Ewing* ' ' ^ and 
by  pressure  data  from  the  Category  1 1/4-scalc  model  tests  reported  by  Schumacher,  et  al.* ' 

For  purposes  of  deriving  the  equations,  the  quasi-static  pressure  is  assumed  to  be  a step  loading 
(zero  rise  time)  to  a constant  value.  If  either  />  or  p is  zero,  the  equations  reduce  to  the  proper 
pressure  or  impulsive  asymptote,  respectively. 

Each  equation  is  based  upon  an  assumed  shape  for  the  final  deformed  state  of  the  struc- 
tural elements.  The  deformation  patterns  used  are  given  in  the  table.  For  some  elements, 
solutions  are  given  for  more  than  one  deformed  shape.  This  is  true,  for  a beam  with  clamped 
ends  which  experiences  bending  deformations  only.  The  first  solution  used  a parabola  as  the 
deformed  shape,  and  the  second  solution  uses  a higher  order  polynomial.  Deformations  pre- 
dicted by  the  first  equation  compare  favorably  to  experiment*'  however,  bending  strains 
in  the  beam  associated  with  a parabolic  deformed  shape  arc  constant.  The  second  equation 
has  net  been  compared  to  expeiiment,  but  the  strain  distributions  produced  by  the  deformed 
shapes  are  more  representative  of  true  beam  behavior.  A third  solution  is  also  given  which  is 
based  upon  the  static  deformed  shape.  Notice  that  the  constants  in  the  equations  derived  for 
the  polynomial  and  for  the  static  deformed  shape  differ  only  slightly. 

Parameters  which  enter  the  equations  in  Table  IV  are  defined  in  Table  V.  Also,  the 
geometry  of  the  element  is  sketched  in  Table  IV  for  additional  clarification.  Because  the 
equations  are  nondimensional  in  their  present  form,  any  consistent  set  of  units  is  permissible. 
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TABLE  rV.  SUMMARY  OF  ENERGY  SOLUTIONS  FOR  ESTIMATING  THE  DEFORMATION  OF 
STRUCTURAL  ELEMENTS  SUBJECTED  TO  UNIFORM  BLAST  LOADING  (Cont’d) 
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I 

I 
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I 
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t 
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TABLE  V.  DEFINITION  OF  SYMBOLS  USED  IN  TABLE  IV 


Symbol 


Definition 


A 

b 

CSBm 

E 

h 

ir 

L 

L,. 

ma 


Mf 

Mp 

N 

P 

Py 

r 

R 

Rl 

Rr 

AR 

w 

X 

X 

y 

Y 

p 

Oy 


beam  cruss-!>cciional  area 
ring  cross-sectional  vea 
loaded  width  of  beam 

circumferential  beam  spacing  in  the  I-beam  cylinder  measured  at  Rp 
material  clastic  modulus 
thickness  of  plate,  or  shell 

specific  reflected  impulse  from  Initial  blast  wave,  plus  reflections 
if  applicable 

length  of  beam  for  which  the  deformation  Is  being  determined; 
length  of  cylinder 

loaded  length  of  the  cylinder  supported  by  a single  ring 

mass  per  unit  area  of  any  additional  material  (non  load  bearing) 
which  is  attached  to  the  sphere  or  dome 

total  mass  supported  by  & ring  (includes  the  ring  mass) 
beam  plastic  moment 

factor  in  tlie  beam  equation;  = 1 for  simple  support,  Af  = 2 for 

clamped  support 

quasi-static  pressure 

axial  yield  force  of  the  beam 

radius  to  arbitrary  point  on  a circular  plate 

mean  radius  of  a sphere  or  cylinder,  radius  of  a circular  plate 

loaded  radius  in  the  cylindrical  shield 

mean  radius  of  the  ring 

radial  expansion  of  the  ring,  dome,  or  cylinder 

lateral  deflection  of  a beam  or  plate  at  point  x,  or  r,  respectively 

center  deflection  of  a beam  or  plate 

distance  along  the  beam  or  plate,  normally  measured  from  the  center 
short  semi-span  of  the  plate 

distance  along  plate  center  line  normally  measured  from  the  plate  center 

long  semi-span  of  the  plate 

material  density 

yield  strength  of  the  material 

Poisson’s  ratio 
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IV.  CATEGORY  I 1/4-SCALE  STRAIN  DATA  ANALYSIS 


McusurciiK-nts  of  strain  in  the  beams,  rings,  foundation,  and  roof  of  the  Category  I 
1/4-si.ale  model  were  made,  by  Schumacher,  et  al.,^'  in  tests  conducted  at  the  BRL. 

Model  dimensions  and  strain  gage  locations  are  shown  in  Figures  5 and  6.  Structural 
details  are  given  in  the  Corps  of  Engineers’  drawing  No.  6003,  “Suppressive  Shield,  Quarter 
Scale  Model,  Category  I.”  One  of  the  principal  objectives  of  the  tests  was  to  determine  the 
structural  adequacy  of  the  shield  and  to  evaluate  the  analytical  methods  used  to  support 
the  shield  design. 

SwRl  was  assigned  the  task  of  analyzing  strain  data  recorded  on  the  beams  and  rings. 
Data  were  provided  by  BRL  in  the  form  of  computer-generated  plots.  The  plots  were  re- 
viewed and  peak  strains  were  extracted  from  records  which  appeared  to  be  consistent  from 
test  to  test  and  with  gages  at  similar  locations.  These  peak  values  were  then  compared 
with  analytical  predictions.  The  comparisons  allowed  us  to  draw  some  conclusions  about 
the  strains  and  analysis  methods  which  are  covered  at  the  end  of  this  chapter. 

A.  Review  and  Summary  of  Experimental  Data 

A brief  description  of  each  test  conducted  on  the  model  is  given  below.  All  charges 
were  spherical  Pentolite,  centrally  located.  Closure  strips  and  liners  referred  to  in  the  test 
descriptions  are  shown  in  Figure  5. 

Test  191:  8.3-lb  charge- 1 /4-scale  model  without  closure  strips  or  liner 

Test  192:  Same  as  Test  191,  but  with  a charge  of  19.3  lb 

Test  193.  19.3-lb  charge-closure  strips  added  to  cover  the  spaces  between  every 

second  pair  of  inside  beams 

Test  194:  19.3-lb  charge  -all  spaces  between  the  inside  beams  covered  with  closure 

strips 

Test  195:  45.7-lb  charge-same  configuration  as  for  Test  194  except  that  additional 

weld  bead  was  added  along  the  sides  of  the  closure  strips,  and  shims  were 
added  to  eliminate  free  travel  by  the  beam  before  contact  with  the  ring 
was  made. 

Test  196:  45.7-lb  charge-weld  repairs  were  made  and  a 24-gage  corrugated  steel 

liner  was  added  inside  the  closure  strips  to  seal  the  shield. 

Test  197:  45.1-lb  charge-ring  repairs  were  made,  and  two  22-gage  corrugated  steel 

liners  were  added  inside  the  closure  strips  to  seal  the  shield. 

We  limited  our  comparisons  between  analysis  and  experiment  to  those  tests  with  all 
closure  strips  installed.  Thus,  comparisons  were  made  with  data  from  Tests  194  through 
197  and  include  comparisons  for  both  the  19.3-lb  and  45.7-lb  charge  weights. 

Figures  7 and  8 are  typical  of  the  strain  plots  received  from  BRL.  Bending  strains  at 
the  base  of  column  1 I 2 (see  Figure  6)  are  given  in  Figure  7 for  Test  192  and  Tests  194 
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FIGURE  7.  BENDING  STRAINS  MEASURED  AT  THE  BASE  OF  COLUMN  112 


FIGURE  8. 


through  197.  The  effect  of  progressive  sealing  of  the  shield  is  apparent  from  the  reduced 
decay  rates.  The  character  of  Test  19S  is  noticeably  different  from  the  others  because  it 
decays  sooner  than  expected  and  has  a negative  value  late  in  time.  Shims  were  inserted  be- 
tween the  beams  and  the  rings  before  this  test,  and  the  change  in  response  character  can  be 
attributed  to  the  change  in  support,  Figure  8 gives  the  history  of  the  principal  strain  com- 
puted from  rosette  No.  85,  which  is  located  at  the  top  of  column  259.  From  these  records 
i*  appears  that  slight  yielding  in  shear  occurred  even  in  the  early  tests.  Reasons  for  the  pro- 
nounced increa.sc  in  shearing  strain  between  Tests  195  and  196  are  not  clear.  Again,  it  may 
have  been  caused  by  firmer  beam  support  after  the  beams  were  shimmed  against  the  rings. 
The  principal  shearing  strain  for  Test  197  (not  shown)  may  have  been  caused  by  firmer  beam 
support  after  the  beams  were  shimmed  against  the  rings.  The  principal  shearing  strain  for 
Test  197  (not  shown)  was  slightly  higher  than  for  Test  196.  A more  detailed  interpretation 
of  the  strain  records  is  given  by  Schumacher,  et  al.^  * ^ ' 

Feak  strains  read  from  the  graphs  of  Figures  7 and  8 and  from  other  graphs  provided  by 
BRL  are  listed  in  Tables  VI,  Vil,  and  VIll.  Table  VI  contains  the  bending  strains  reduced 
from  gage  pairs  on  the  inside  and  outside  flanges  of  the  beams.  Location  of  the  gages  above 
the  top  surface  of  the  foundation  is  included  in  the  table,  and  the  placement  of  the  gages  on 
the  beam  is  shown  in  Figure  6.  Principal  shearing  strains  in  the  beams  are  given  in  Table  VII 
and  peak  tensile  strains  in  the  rings  are  given  in  Table  VIII.  In  all  cases,  the  strain  data  for 
Test  191  (the  8.3-lb  test)  and  for  Test  193  (the  19.3-lb  test  with  only  half  of  the  closure 
strips  installed)  have  been  omitted.  Test  192  was  included  to  indicate  the  effect  which  the 
closure  strips  have  on  the  peak  strains.  No  consistent  difference  is  apparent  in  the  beam 
strains.  Ring  strains  are  increased  slightly. 

TABLE  VI.  MAXIMUM  ."ENDING  STRAINS*  OBTAINED  FROM  BRL  RECX)'.<DS 
OF  CATEGORY  I 1/4  SCALE  MODEL  TESTS 


Locatlont 

(In.) 

Gage 

Nos. 

19.31b  Charge 

45.7-lb  Chai,>t; 

Test  192 

Test  194 

Test  195 

Test  196 

Test  197 

6-3/8 

50-51 

1050 

840 

700 

llf'U 

1250 

56-57 

1175 

1200 

- 

8/5 

- 

66-67 

- 

- 

- 

— 

- 

72-73 

- 

- 

- 

- 

- 

34-1/2 

58-59 

-975 

-1000 

-1450 

- 

- 

45 

60-61 

-1850 

— 

— 

— 

74-75 

_ 

-1350 

-1165 

— ■ 

-1600 

68-69 

- 

- 

- 

- . 

- 

55-1/2 

62-63 

-1000 

-1000 

-900 

- 

-1450 

113-.5/8 

54-55 

1300 

1225 

1090 

1155 

1135 

64-65 

1100 

1240 

- 

11  sot 

- 

70-71 

- 

- 

- 

— ■ 

— 

76-77 

1350 

1450 

1150 

— 

*Stri>lni  are  in  >iin/ln.i  (+)  tepreiemi  teniion  on  iniide  vds.'  of  the  beam. 
tMnanired  /roni  top  of  floor  slab  (lee  Figure  6). 

^Eitimatrd  from  one  gage  only. 
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TABU-  VII.  MAXIMUM  SHEAR  STRAIN*  OBTAINED  FROM  BRL  RECORDS 
OF  CATEGORY  I 1/4  SCALE  MODEL  TESTS 


Location^ 

(in.) 

Rosette 

No. 

19.3-lb  Charge 

45.7-lb  Charge 

Test  192 

Test  194 

Test  195 

Test  196 

Test  19'? 

34-1/2 

78 

1700 

2650 

3250 

81 

- 

- 

- 

- 

- 

113-S/8 

80 

- 

- 

- 

83 

21  SO 

2000 

2500 

7500 

- 

84 

-- 

1550 

2400 

- 

85 

2800 

2350 

7000 

14000 

14810 

‘Slruins  are  In  uln/ln 

t Measured  from  top  of  flnot  slab  (see  Figure  6). 


TABLE  Vlll.  MAXIMUM  RING  STRAINS*  OBTAINED  FROM  BRL  RECORDS 
OF  CATEGORY  I 1 /4-SCALE  MODEL  TESTS 


Locationt 

Gage 

No. 

19.3- 

lb  Charge 

45.7-lb  Charge  | 

Test  192 

Test  194 

Test  195 

Test  196 

Test  197 

Tup  Ring 

90 

1060 

2580 

2400 

2200 

91 

-1430 

3100 

2400 

4200 

Middle 

92 

1450 

1600 

_ 

3300 

3400 

Ring 

93 

1150 

2300 

2800 

2100 

3100 

'Strains  arc  in  pln/ln. 
tScB  Figure  6. 


B Comparisons  with  Analytical  Predictions 

For  comparison  with  the  measured  strains  of  Tables  VI,  VII,  and  Vlll,  strains  in  the 
beams  and  rings  were  computed  usin^!  both  the  approximate  energy  procedures  described 
in  Section  II  and  finite-element  methods.  Application  of  the  energy  methods  to  compute 
structural  response  has  been  reported  in  References  1,  4,  7,  9,  and  10,  and  in  earlier  chapters 
of  this  report,  so  these  energy  methods  will  be  used  here  with  a minimum  of  explanation. 
Application  of  the  finite  element  program  to  compute  structural  response  will  be  described 
more  thoroughly. 

/.  Predictions  Using  Approximate  Energy  Methods 
a.  Uncoupled  Solution 

The  blast  loads  associated  with  a confined  explosion  as  in  a suppiessive  struc- 
ture are  not  the  same  as  the  blast  loads  associated  with  -an  unconfined  explosion.  Initially  in 
a confined  explosion,  a shocK  wave  is  propagated  out  away  from  a source;  however,  because 
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of  the  walls  in  the  container,  this  initial  shock  is  reflected  many  times  until  through  various 
dissipation  mechanisms,  the  air  is  heated  and  a static  pressure  buildup  of  very  long  duration 
results.  In  the  suppressive  structures  program,  this  multiple  loading  mechanism  was  mathe- 
matically modeled  by  treating  the  initial  shock  wave  as  a delta  function  (as  an  impulse)  and 
the  long  duration  buildup  of  internal  pressure  as  a constant  static  pressure. 

The  generalized  solution  for  the  response  of  a beam  in  bending  to  this  type  of 
loading  is  given  in  Figure  4 of  Chapter  II.  As  in  all  solutions  for  structural  elements;  ideal 
boundary  conditions  are  assumed,  and  coupling  between  the  element  and  its  supporting 
structure  are  neglected.  For  the  Category  1 model,  the  assumption  of  clamped-clamped 
boundaries  for  the  beams  is  most  appropriate. 

Beams  in  the  model  are  S3  X 5.7.  Properties  of  the  beam  cross-section  arc 
obtained  from  the  Steel  Construction  Handbook  and  the  remainder  of  the  beam  geometry 
from  Figure  5.  Loads  on  the  beam  are  produced  by  the  initial  reflected  blast  wave  and  by 
the  subsequent  quari-static  pressure.  Impulse  in  the  blast  wave  is  obtained  from  the  data 
presented  by  Baker.* ' The  quasi-static  pressure  is  determined  using  the  procedure  developed 
by  Baker,  et  al.* ' ' Compiling  the  input  data  for  Figure  4 from  the  above  sources,  we  have 

Geometry: 

A = 1.67  in.^  (area  of  the  beam) 

b=  1.43  in.  (loaded  width-less  than  flange  width  because  of  overlap) 

H = 3.0  in.  (beam  depth) 

z~  1.95  in.**  (beam  plastic  section  moduhis) 

/=  1.65  in.^  (beam  elastic  section  modulus) 

Oy  = 45,000  psi  (material  yield  strength) 

£ = 30  X 10‘  psi  (elastic  modulus) 
p = 7.33  X 10'“  lb-sec^ /in. '*  (material  density) 

Loading  Parameters: 

19.3-lb  Charge  45.7-lb  Charge 

Reflected  impulse,  L = 0.20  psi-sec  0.48  psi-sec 

Quasi-static  pressure,  Pqj  = 73  psi  190  psi 

A yield  strength  greater  than  the  minimum  guaranteed  for  the  material  was  chosen  because 
limited  measurement  on  structural  steel  conducted  during  the  suppressive  structures  program 
showed  yield  values  between  40.000  psi  and  50,000  psi. 

Maximum  ueam  strain  for  the  45.7-lb  charge  weight  was  computed  in  Chapter 
II,  Part  D,  as  an  example  for  the  use  of  Figure  4.  The  calculated  strains  given  in  Table  IX  were 
obtained  using  the  same  procedure  for  the  loading  produced  by  the  19.3-lb  charge. 
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TABLE  IX.  COMPARISON  OF  PEAK  BEAM  BENDING  STRAINS: 
EXPERIMENT  TO  UNCOUPLED  ENERGY  SOLUTIONS 


Charge  Weight 

Calculated 

Measured 

(lb) 

m 

(|je) 

19.3 

3,730 

1450 

45.7 

12,700 

1600 

For  the  19.3-lb  charge,  data  from  Test  194  rather  than  Test  192  were  used.  The  peak  strain 
occurred  at  the  top  of  an  inside  beam.  For  the  45. 7-lb  charge,  the  maximum  measured  strain 
was  midway  between  two  supports  (two  rings)  on  the  same  beam.  The  gages  which  recorded 
the  maximum  strain  for  the  19.3-lb  charge  were  no  longer  operative  for  the  tests  with  the 
45.7-lb  charge,  so  bending  strains  at  the  top  of  the  beam  might  have  been  greater.  Even  if 
this  were  so,  it  is  apparent  that  the  predicted  strains  are  much  larger  than  measured  values. 
This  is  particularly  true  for  the  45.7-lb  charge. 

The  predictions  of  Table  IX  were  obtained  by,  in  effect,  treating  the  rings  as 
rigid  and  absorbing  al-  of  the  energy  produced  by  the  blast  in  the  beams  only.  Obviously, 
these  predictions  ate  much  too  conservative  to  be  of  practical  use  in  shield  design.  Similarly, 
ring  response,  computed  for  the  assumption  of  rigid  beams,  is  much  too  great.  Thus,  the 
energy  absorbed  must  be  properly  divided  between  the  beams  and  the  rings,  and  this  requires 
a coupled  solution. 


b.  Coupled  Solution 

A solution  for  the  coupled  response  of  two  structural  components,  one  sup- 
ported by  the  other,  is  described  in  Chapter  II  and  Reference  10.  It  assumes  rigid-plastic 
behavior  with  loads  applied  to  the  supported  structure  only.  The  loading  is  comprised  of 
an  initial  impulse  and  a quasi-static  pressure  as  for  the  uncoupled  solution  described  previously. 
For  beams  supported  by  rings  (hoops)  as  in  the  Category  I .shield,  the  response  of  the  system 
is  given  by  Eq.  (64)  of  Chapter  II.  This  equation,  repeated  here  for  convenience,  gives  the 
motion  of  the  center  of  the  beam,  Wo,  relative  to  the  rings. 


, pbi^  , NphAhU  \o/i/4a8  / „ . OhPbAbi 

— 1,0 -F 1 i = nTT  T 

PbAiiObZWo  NobZ  4itpi,AhH  pri  \ ObPhzR 

V 2ai,Ah/ 


flNabZ 
OhPbAbi^  \ohAhi 


/2A/o^  _ j \ 

\ohAH«.  / 
V 2otA,/ 


Once  Wo  is  found,  Eq.  (59),  also  taken  from  Chapter  11,  gives  the  total  radial 
motion,  AR,  of  the  ring  in  terms  of  Wo. 


w„  ^ \Npi,Abi/  \ 2of,Ah/ 

l\  a„Ahi/  J 
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Based  on  the  maximum  deflections,  Wg  and  AR  of  the  beam  and  ring  respectively,  the  residual 
strains  aie  computed  from  Eqs.  (65u)  and  (65b). 


(65a) 


et 


it^WoH 


(65b) 


To  compute  residual  strains  in  the  Category  I 1 /4-scale  model  using  these 
equations,  the  following  parameters  are  required: 


Ob  - Oh  = 45,000  psi  (material  yield  strength) 

Pb  ~ Ph  - ‘^•33  X 10‘*  lb-sec^ /in. ^ (material  density) 

Ab  = 1.67  in.^  (beam  cross-section  area) 

2 = 1.95  in.^  (beam  plastic  section  modulus) 

N = 296  (no.  of  beams  supported  by  the  rings) 

6 = 1.43  in.  (beam  loaded  width- less  than  fiange  width  because  of  overlap) 
8 = 30  in.  (beam  length) 

H - 3.0  in.  (beam  depth) 

Ah  = 5.625  in.^  (one-half  of  the  ring  area) 

/?  = 72.25  in.  (radius  to  center  of  the  ring) 
r = 67.25  in.  (internal  radius  of  the  shield) 


Blast  loads  are  the  same  as  listed  on  page  62. 


Results  obtained  with  the  coupled  solution  are  compared  to  the  experiment 
in  Table  X.  Both  beam  and  ring  strains  are  included.  Calculated  strains  are  residual,  whereas 
the  measured  strains,  are  peak.  It  is  not  appropriate  to  simply  add  the  elastic  strain  at  yield 
(Cy  = Oy/E  - 1550  Min. /in.)  to  the  residual  strain,  but  the  peak  strain  would  lie  between  the 
residual  strain,  Cr,  and  the  residual  plus  elastic  strain  (6;^  + 6c).  With  this  in  mind,  the  agree- 
ment between  the  calculated  and  measured  strains  is  good.  Very  little  yielding,  if  any, 
occurred  in  the  beams,  and  this  is  predicted  by  the  analysis.  Predicted  ring  strains  are  higher 
than  measured  for  the  45.7-lb  charge.  One  reason  for  this  may  be  that  the  nine  Individual 
bands  from  which  the  ring  was  assembled  were  not  in  good  contact.  This  means  that  higher 
strains  occurred  in  the  inner  bands  than  in  the  outer  bands  where  the  measurements  were 
made.  Also,  many  rings  broke  near  the  welds,  which  may  indicate  an  embrittlement  problem. 
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TABLE  X.  COMPARISON  OP  BEAM  AND  RING  STRAINS: 
EXPERIMENT  TO  COUPLED  ENERGY  SOLUTIONS 


19.3-lb  Charge 

45.7.1b  Charge 

Structural 

Calculated 

Measured 

Calculated 

Measured 

Component 

Residual 

Uie) 

Residual 

(pe) 

(fie) 

Beams 

200 

1450 

390 

1600 

Rings 

97S 

2300 

6.880 

4200 

As  can  be  seen  by  comparing  the  results  in  Tables  IX  and  X,  the  use  of  coupling 
is  to  greatly  reduce  the  calculated  estimates.  Although  the  estimates  are  still  not  in  perfect 
agreement,  the  coupled  predictions  are  much  closer  to  measured  values  than  are  the  uncoupled 
predictions.  Even  better  prediction  should  be  achieved  with  finite  element  methods. 

2 Predictions  by  Finite  Element  Methods 

In  addition  to  the  energy  solutions,  a small  two-dimensional  (2-D)  finite  element 
computer  program,  developed  specifically  for  the  dynamic-transient  analysis  of  beams  to 
blast  loading,  was  used  to  predict  structural  response  of  the  CategOiy  1 1 /4-scale  model.  The 
equations  of  motion  ars  cast  in  finite  difference  form  and  integrated  step-wi<.e  in  time  using 
a predictor-corrector  method.  The  program  treats  material  behavior  as  bilinear  with  hysteretic 
recovery,  including  the  effects  of  strain  rate  and  strain  hardening.  Calculations  of  the  response 
were  made  with  and  without  strain  rate  to  show  the  significant  effect  strain  rate  has  on  dis- 
placements and  strains.  Strain  hardening  effects  were  found  to  be  negiigible  and  are  not  dis- 
cussed further. 

The  finite  element  model  used  for  the  calculations  is  shown  in  Figure  9.  It  is  drawn 
schematically  alongside  one  of  the  beams  from  the  cylinder.  The  structure  is  idealized  by 
assuming  axial  symmetry  and  symmetry  about  a horizontal  plane  through  the  center  of  the 
shield.  This  amounts  to  neglecting  the  presence  of  the  door  on  the  response  of  the  beams 
away  from  its  immediate  vicinity,  and  neglecting  axial  loading  and  response.  Although  there 
is  some  connection  between  the  beams  through  the  closure  strips,  which  are  tack  welded  to 
the  face  of  the  inner  beam  and  to  the  studs  on  the  outer  beam  (see  Figure  5),  the  beams  are 
assumed  to  respond  individually.  Thus,  only  the  lower  half  of  one  beam  is  represented  in 
the  model. 

Bilinear  springs  have  been  used  to  model  the  rings  and  radial  flexibility  in  the 
foundation.  The  resistance  of  the  rings  is  divided  into  two  parts  to  spread  the  support  over 
the  ring  width.  Guided  end  conditions  permit  lateral  displacements,  but  no  rotation,  and 
the  first  beam  element  above  the  foundation  has  been  sized  to  represent  the  combined  bending 
stiffness  of  the  slab  foundation  and  the  bolted  connection.  Although  the  closure  strips  may 
contribute  to  the  beam  stiffness,  they  are  not  welded  continuously  to  the  beams,  and  so  this 
contribution  is  ignored;  however,  mass  of  the  closure  strips  and  the  liners  is  added  to  the 
beam  mass. 
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HGURE  9.  SCHEMATIC  OF  THE  FINITE  ELEMENT 
MODEL  OF  THE  STRUCTURE 
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As  for  the  energy  solutions,  the  pressure  loading  was  assumed  to  be  uniformly  dis- 
tributed over  the  beam.  The  loading  is  composed  of  an  initial  Impulse  and  near  constant 
(quasi-static)  pressure.  Two  load  cases  were  analyzed.  They  are  the  same  as  listed  on  page 

19,3-ib  charge:  i,  = 0.20  psi-sec 

P</i  = 73  psi 

45.7-lb  charge:  1,  ~ 0.48  psi-sec 

Pqs^  190  psi 

Basic  beam  and  ring  parameters  used  in  computing  input  for  the  program  are: 

• S3  X 5.7  beam  geometry 

• mass  per  unit  length  (including  closure  strips  and  liners), 

A = 1.46  X 10'^  Ib-sec/in.^ 

• foundation  effective  modulus,  tc  = 4.5  X 10‘  p.si 

• beam  plastic  moment,  = 87.750  in. /lb 

• beam  elastic  modulus,  £■  = 29  X 1 0®  psi 

• beam  plastic  modulus  = l%of  elastic  modulus 

• beam  loaded  width,  i = 1 .43  in. 

Although  the  details  will  not  be  given,  stiffness  of  the  reinforcing  gussets  at  the  foundation 
and  of  the  spacers  at  the  rings  was  added  to  the  beam.  Also,  estimates  were  made  of  the 
bending  stiffness  of  the  beam  attachment  to  the  foundation.  Strain  rate  effects  are  based 
on  the  well  known  relationship^ for  mild  steel  given  in  Eq.  (153) 

o<,  = a,  [>+(4-)°']  (153) 

where  oj  is  the  dynamic  yield  stress,  and  e is  the  strain  rate. 


Results  for  the  19.3-lb  charge  are  given  in  Figure  10.  Strain  rate  effects  are  included. 
This  figure  shows  the  distribution  of  maximum  values  of  displacement,  bending  strain,  and 
shear  stress  along  the  beam.  The  shear  stress  is  based  on  the  shear  reactions  in  the  beam  rather 
than  on  shearing  strains,  and  is  computed  as 


(154) 


where  V is  the  shear  reaction  in  the  beam,  d is  the  beam  depth,  and  is  the  beam  web  thick- 
ness. Data  points  are  shown  as  circles.  Three  observations  can  be  made  from  this  comparison: 


( I ) The  distribution  of  peak  .strains  in  the  beam  is  irregular  and  not  what  one 
would  expect  if  only  a section  of  the  beam  between  two  rings  were  analyzed. 
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Radial  Displacement,  6^.  (in.  ) 


FIGURE  10.  DISTRIBUTION  OF  MAXIMUM  DEFLECTION.  ^.TRAIN 
AND  SHEAR  STRESS  VALUES  (19.3  lb  chaiirl 


(2)  Computed  results  for  bending  strains  underestimate  the  measured  values  for 
this  charge  weight  (at  the  measurement  location). 

(3)  If  the  computed  distribution  of  strain  is  accurate,  the  measurement  point 
may  not  have  recorded  the  maximum  values. 

Figure  1 1 gives  similar  results  for  the  45.7-lb  charge,  The  same  observations  can 
be  made  for  this  case  as  were  made  for  the  smaller  charge  excepi  that,  in  this  instance,  the 
computed  strains  are  higher  than  measured  values.  Also,  for  this  charge  weight  the  shearing 
stresses  exceed  the  shear  yield  stress  of  the  material  at  many  points  along  the  beam.  We 
estimate  that  the  allowable  shear  stress  values  for  the  material  are: 

a,y  = 21,600-27,000  psi 

Oju  = 34,800—49,800  psi 

These  compare  to  peak  calculated  value:;  of  44,000  psi,  which  indicates  that  substantial 
yielding  of  the  beam  in  shear  should  have  occurred,  and  perhaps  the  beam  was  near  shear 
failure. 


Although  strain  comparisons  cannot  be  made  because  shearing  strains  were  not 
computed  by  the  program,  Table  XI  lists  measured  shear  strains  versus  the  computed  stress 
The  measured  results  confirm  that  substantial  yielding  occurred  in  shear  for  the  45.7-lb 
charge.  Shear  strains  slightly  beyond  yield  were  also  recorded  for  the  19.3-lb  test.  These 
results,  plus  the  results  in  Figure  1 1 , indicate  that  the  beam  is  not  properly  proportioned 
for  1 ,;;hly  impulsive  short  duration  loading.  Higher  shear  area  relative  to  bending  stiffness 
is  needed  to  bring  the  shearing  and  bending  strains  into  better  balance.  This  observation  is 
true  for  alrr.ost  all  commercially  available  beams  b.-.cause  they  have  been  designed  for  high 
ef  ciency  under  static  loads. 


TABLE  XI.  BEAM  SHEAR  DATA 


19.3-lb 

45.7-16 

Location 

Meas. 
Strain 
Test  194 

Calc. 

Stress 

F.E.  Prog 

Measured  Strain 

Calc. 

Stress 

F.E.  Prog 

— 

Test  195 

Test  )96 

Test  i97 

Outer  beam  34-1/2" 
from  foundation 

2,650* 

18,0001 

3,250 

-- 

-- 

35,000 

Outer  beam  6-1/4" 
from  roof 

1.550 

20,000 

" 

- 

44,000 

Inner  beam  6-1/4" 
from  roof 

2,000 

20,000 

- 

44,000 

Inner  beam  6-1/4" 
from  roof 

2,350 

20,000 

14,800 

44,00(' 



■Strain  M In./in. 
(Stress  psi. 
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FiCURE  11.  DISTRIBIITICN  OF  MAXIMUM  DEFLECTIOM,  STRAIN 
AND  SHEAR  STRESS  VALUES  (4S.7  lb  chaifo) 
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Only  results  which  included  strain  rate  effects  were  presented  In  Figures  10  and  1 1. 

The  effects  of  strain  rate  on  the  overall  structural  response  are  best  demonstrated  by  the  variation 
in  ring  strains.  Measured  and  calculated  strains  are  given  in  Table  XII.  Calculated  values  include 
those  obtained  by  energy  formulas  as  well  as  by  the  finite  element  program.  For  the  energy 
solutions,  both  coupled  and  uncoupled  results  are  included.  The  coupled  solution  was  obtained 
from  Eqs.  (59)  and  (64);  the  equation  for  uncoupled  response  is  given  in  Table  IV  of  Chapter  III. 

TABLE  XII.  MEASURED  AND  CALCULATED  STRAINS  IN  THE  RINGS 


19.31b 

45.7  lb 

Test  194 

Test  195 

Test  196 

Test  197 

Top  Ping 

Measured 

Gage  90* 

Gage  91 

11001 

2200 

3100 

2400 

2400 

2200 

4200 

F.E.J  Solution 

No  S.R.** 

With  S.R. 

5315 

1840 

13630 

5110 

Energy  Formulas 

No  coupling 

With  coupling 

3100 

2460 

18490 

6320 

Center  Ring 

Measured 

Gage  92 

Gage  93 

1600 

2300 

“ 

3300 

2100 

3400 

3100 

F.E.  Solution 

No  S.R. 

With  S.R. 

10360 

3000 

21530 

9740 

Energy  Formulas 

No  coupling 

With  coupling 

3100 

2490 

18490 

6824 

'Refer  to  Figu>  t 7 for  gage  locations, 
f Strains  ate  in  ;i  itt./in. 

^Finite  element  2-D  program. 

"Strain  rate  effects. 

The  most  obvious  conclusion  to  be  reached  from  the  results  in  Table  XII  is  that 
the  calculated  values  overestimate  measured  ring  response  in  every  case.  For  the  two  com- 
parisons available,  it  appears  that  the  difference  between  computed  and  measured  values 
increases  as  charge  weight  increases:  hence,  the  analytical  results  should  always  be  conservative. 
The  effects  of  strain  rate  and  coupling  are  pronounced,  reducing  calculated  values  by  approxi- 
mately a factor  of  three  for  the  large  charge  weight.  Still,  the  calculated  values  are  far  above 
measured  values.  It  should  be  noted  that  the  energy  solution  does  not  include  strain  rate  or 
strain  hardening  effects  and  yet  gives  a result  which  is  the  same  order  of  magnitude  as  the 


71 


more  rigorous  finite  element  method.  One  reason  for  this  is  that  the  energy  solution  forces 
plastic  deformation  in  the  beams,  whereas  the  finite  element  solution  does  not.  It  is  apparent 
from  Figure  I 1 that  predicted  plastic  straining  in  the  beams  was  minimal  so  that  energy  absorp- 
tion by  the  beams  was  small. 

Reasons  for  the  disagreement  between  predictions  and  measurements  of  ring 
response  are  not  deal . However,  two  factors  may  have  affected  the  results; 

(1 J The  closure  strips  shown  in  Figure  5 were  tack  welded  to  the  inner  beams  and 
provided  some  circumferential  restraint.  Thus,  the  closure  strips  acted  in  a 
manner  similar  to  rings,  but  were  distributed  over  the  full  height  of  the  shield. 
This  effect  was  neglected  in  the  response  calculations. 

(2)  Properties  of  the  rings  material  were  never  measured.  U proved  to  be  a brittle 
material  because  fractures  were  noted  during  the  testing.  The  calculatior.s 
were  performed  assuming  a mild  steel,  which  may  not  represent  true  ring 
properties. 

One  final  comparison  of  measured  and  calculated  results  is  given  in  Table  XHI.  It 
compares  peak  measi'*'ed  strains  in  the  beam  to  strains  and  shearing  stresses  computed  by 
both  the  energy  and  tinite  element  methods.  From  this  comparison,  it  is  obvious  that  both 
the  finite  element  solution  with  strain  rate  effects  and  the  coupled  energy  solution  give  results 
which  are  consistent  with  the  measured  values.  Additionally,  the  finite  element  method  pro- 
vides useful  information  on  shear  in  beams  by  providing  shear  reactions  which  predict  the 
presence  of  large  shearing  strains.  From  the  finite  element  solution,  it  is  also  apparent  that 
strain  rate  effects  play  little  part  in  the  beam  bending  response.  By  stiffening  tne  ring,  how- 
ever, strain  rate  does  cause  higher  shear  stresses.  This  stiffening  of  the  rinj,  was  demonstrated 
in  Table  XII  by  the  sharp  reduction  in  the  ring  strain. 


TABLE  XJII  COMPARISON  01'  PEAK  VALUES  IN  THE  BEAMS 


19.3  lb  charge 

45.7  lb  charge 

Bending 

Strain 

H in. /in 

.Shear 
Strain 
(i  in  ./in. 

Shear 

Stress 

psl 

Bending 

Strain 

M in. /in. 

Shear 

Strain 

M in. /in. 

Shear 

Stress 

psi 

Measured 

1450 

2650 

1600 

14,800 

F.E.  Method 

No  SR 

With  SR 

1400 

20,000 

3010 

3000 

36,900 

44,000 

Energy  Solution 
Uncoupled 
Coupled 
(Residual) 

3730 

200 

12700 

390 

C.  Conclusions  from  the  Strain  Data  Analysis 

The  conclusions  drawn  from  these  comparisons  of  analysis  and  experiment  are: 

( 1 ) Strain  data,  provided  by  BRl.  in  the  form  of  computer  generated  plots,  appear  to 
be  good.  The  data  are  surprisingly  consistent  from  test  to  test  for  a structure  which 
was  being  altered  continuously  and  undergoing  plastic  strains.  Exceptions  occur, 
and  these  can  usually  be  explained  by  major  alterations  to  the  structure,  such  as 
shintming  the  beams  prior  to  Test  195. 

(2)  Measured  beam  bending  strains  are  lower  than  expected  and  did  not  increase  sig- 
nificantly with  an  increase  in  charge  weight. 

(31  Shear  strains  in  the  beam  web  were  high.  Yielding  occurred  as  predicted  by  high  . 
shear  forces  which  were  calculated  for  the  beams. 

(4)  Ring  strains  were  much  lower  than  predicted,  particularly  for  the  large  charge 
weight. 

(5)  The  comparisons  show  that  the  energy  solution  for  coupled  response  can  give 
reasonable  answers,  which,  if  in  error,  tend  to  overestim.ate  the  strains. 

(6)  The  finite  element  computer  solutions  show  that  strain  rare  should  not  be  ignored 
in  structural  componenis  which  respond  primarily  in  a stretching  mode,  such  as  for 
the  rings.  If  these  effects  are  ignored,  the  solution  is  an  upper  bound,  but  it  is  too 
conservative  to  be  useful  if  reasonably  accurate  predictions  are  desired. 


V.  PRESSURES  FROM  BURNING  PROPELLANT 
IN  VENTED  CHAMBERS 

The  thrust  of  this  effort  was  to  determine  theoretically  the  peak  quasi-static  pressure 
produced  by  burning  M 10  propellant  in  closed  and  vented  chambers.  Combustion  equations 
were  written  to  define  the  heat  generation  and  products  of  the  burning  propellant  which  act 
to  increase  the  pressure  in  the  shield.  Temperature  and  pressure  increases  in  the  shield  pro- 
duced by  burning  the  propellant  are  offset  by  radiatipn  to  the  shield  walls  and  gas  flow  from 
the  sfiield  through  the  vented  panels.  The  equations  describing  these  processes  were  coded 
into  a computer  program  called  GASSS  which  performed  the  calculations  for  M 10  propellant 
burned  in  a Category  V shield. 

A.  Combustion  Equations 

M 10  propellant  is  98'/c  nitrocellulose  with  13.15%  nitration.  Other  m^or  constituents 
are  potassium  sulfate  and  diphenylamine.  Since  nitrocellulose  is  the  most  active  element  in 
this  propellant  and  since  most  of  the  propellant  is  nitrocellulose,  the  program  assumes  that 
MIO  is  all  nitrocellulose.  Thus,  our  calculations  for  energy  produced  from  burning  the  pro- 
pellant will  be  an  upper  limit.  Based  on  13.15%  nitration  of  nitrocellulose  and  information 
in  the  literature  on  its  products  of  combustion,^ the  reaction  equations  for  burning  nitro- 
cellulose in  the  absence  of  air  are: 
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05-*1.556  C0,  + 4.407  CO  + 1.001  Hj  + 1.29.7  N,  +2.622  HjO 


+ 265.0  kilocalories  (155) 

Secondary  reactions  involving  the  combustion  of  CO  and  Hj  in  air  are; 

CO  + 1 /20,  -CO,  + 67 , 63  kilocalories  (156) 

and 

H,  + l/20,-H,0  + 57.80  kilocalories  (157) 

The  program  was  written  so  that  during  each  time  step,  the  products  and  energy  of  the 
reaction  are  calculated  first  without  the  presence  of  oxygen.  Then,  a calculation  is  performed 
to  determine  the  amount  of  oxygen  available.  If  there  is  sufficient  oxygen  present,  the  si:con- 
dary  reactions  take  place.  When  the  oi.ygen  in  the  chamber  is  depleted,  only  the  primary 
reaction  takes  place.  In  vented  and  closed  chambers,  it  is  assumed  that  no  oxygen  enters 
from  outside  the  structure  due  to  positive  pressure  inside  the  chamber.  The  computer  pro- 
gram also  records  the  quantity  of  different  gases  present  in  the  chamber  at  all  times. 

B.  Gas  Flow  Equations 

Determination  of  the  pressure  history  for  venting  in  a suppressive  structure  was  accom- 
plished by  assuming  quasi-steady,  isentropic  flow  of  a perfect  gas.  A check  on  the  reliability 
of  the  perfect  gas  assumption  was  made  by  assuming  complete  combustion  in  an  unvented 
chamber  with  no  energy  loss,  and  thermal  equilibrium  in  the  gases.  Pressure  was  calculated 
using  the  ideal  gas  law  and  also  using  Van  der  Waal  s equations  for  non-ideal  gases,  and  the 
results  were  compared.  For  the  pressures  and  molecular  quantities  involved,  the  ideal  gas 
law  was  found  to  be  quite  accurate. 

For  the  calculated  pressure  dissipation,  the  suppressive  structure  was  idealized  as  two 
compartments  separated  by  some  vent  area.  If  the  ratio  of  atmospheric  pressure  to  chamber 
pressure  is  less  thaii 


7+1/ 

then  supers-j.nic  flow  equations  are  used.  Otherwise,  subsonic  flow  equations  are  utilized. 
The  program  causes  the  gas  to  vent  proportionally  to  the  quantity  of  each  gas  present. 

If  the  flow  is  supersonic,  the  mass  flow  rate  is  determined  by^ 
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where 

A - vent  area  of  passage  out  of  chamber 
Co  - discharge  coefficient  of  area  A 
Pi  = compartment  pressure 
R = universal  gas  constant 
T - gas  temperature 
g = gravity  constant 
7 = ratio  of  specific  heats 

If  the  flow  is  subsonic,  then 

At 

where 

Pi  = density  of  gas  in  the  chamber 
p2  = pressure  outside  the  chamber 

The  rate  of  temperature  change  is  given  by 
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where 

m = mass  of  the  gas  in  the  compartment,  and  tlie  rate  of  pressure  change  is  given  by 
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where 


V - volume  of  the  compartment. 
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C.  Effects  of  Burning  Rates  and  Radiant  Heat  Loss 

Insufficient  data  were  available  to  determine  accurately  the  burning  rates  of  piles  of  MIO 
propellant.  A telephone  conversation  with  Mr.  William  Seals/* of  Picatinny  Arsenal,  pro- 
duced the  following  results  for  MIO: 


Burning  Rate  (in./sec) 

Pressure  (psi)  Single  Perforation  Multi-perforation 

5.000  1.49  1.38 

9.000  2.24  2.07 

11,000  2.58  2.18 

These  burning  rates  are  obviously  appropriate  for  interior  gun  ballistics,  but  are  not  appropriate 
for  our  purposes  because  the  pressures  in  the  suppressive  structure  for  the  worst  case  (no  energy 
loss  and  instantaneous  burning)  are  much  smaller  than  the  pressures  in  the  table  above.  Measure- 
ments that  Dr.  Gary  McKown^*^*  had  taken  when  burning  large  quantities  (500-1000  lb)  of 
smoke  mixes  indicated  burning  rates  of  about  10  Ib/sec.  Measurements  which  he  had  taken 
while  burning  10-  and  50-lb  quantities  of  magnesium  and  sodium  nitrate  illuminant  produced 
burning  rates  varying  from  4 to  100  Ib/sec,  depending  on  the  quantity  burned  and  the  degree 
of  confinement. 

In  view  of  this  uncertainty  in  the  burning  rate,  various  rates  were  used  in  the  calculations 
to  show  their  influence  on  peak  pressures  in  the  Category  V shield.  Initially,  burning  rate  was 
defined  in  units  of  ib/sec,  and  the  pressure  histories  in  the  shield  were  computed  using  various 
quantities  of  propeliant,  several  different  burning  rates,  and  various  effective  vent  areas.  The 
resulting  pressures  seemed  to  be  unrealistically  high,  and  this  approach  appeared  to  be  inaccurate 
since  the  burning  rate  should  depend  on  surface  area  and  ambient  pressure,  and  there  should  be 
radiant  energy  loss  depending  on  the  buniing  surface  temperature  and  area.  Thus,  the  computer 
program  was  modified  to  incorporate  a linear  burning  rate  and  radiant  heat  lo.ss.  Because  of  the 
relatively  low  peak  pressures,  we  neglected  pressure  dependency  of  the  burning  rate. 

To  convert  a linear  burning  rate  to  a mass  burning  rate  and  burning  surface  area  required 
that  a burning  geometry  be  assumed.  Further,  it  was  necessary  that  the  burning  surface  be 
small  initially  and  increase  with  burning  time,  which  is  characteristic  of  point  ignition.  The 
simplest  geometry  which  could  be  assumed  was  a sphere  witli  point  ignition  at  the  center.  To 
permit  radiation  the  sphere  was  treated  as  two  hemispheres.  For  this  geometry  the  effective 
increase  in  surface  area  with  time  is  equal  to  that  of  a radially  expanding  sphere.  The  largest 
burning  surface  obtained  in  using  this  procedure  is  nearly  identical  to  that  of  a cube  of  pro- 
pellant of  equ.al  weight  with  one  face  unexposed;  that  is,  on  the  floor.  For  example,  equating 
volumes  of  the  sphere  and  cube: 


8^ 


(163) 


where  R = radius  of  sphere  (maximum) 
fi  = edge  of  cube 
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Solving  Eq.  (163)  for  /?,  one  has 


t 


• I 


The  maximum  surface  area  of  the  sphere  is 

4irr^  = 47r  ( — 

\47r 

which  is  nearly  the  area  of  the  five  exposed  sides  of  a cube  of  identical  weight. 

The  above  assumptions  on  bulk  propellant  burning  and  increase  in  burning  surface  were 
made  to  simplify  a very  complex  burning  process  and  provide  for  a representative  increase  in 
burning  surface  area  and  attendant  increasing  radiation  losses  after  ignition.  The  use  of  a 
spherical  geometry  to  describe  the  burning  process  offers  the  advantage  of  simplicity.  !t  also 
allows  for  point  ignition  and  growth  in  burning  surface  to  a final  value  which  is  nearly  the 
same  in  area  as  the  exposed  surfaces  of  a cube  of  equivalent  weight.  Thus,  the  simplified 
spherical  geometry  very  nearly  approximates  the  characteristics  of  the  rrctual  burning  process. 

The  computer  program  was  changed  to  incorporate  a radial  burning  rate  (uiiits  of  in. /sec) 
for  the  bulk  propellant.  By  assuming  negligible  pressure  dependency  (relatively  low  peak 
pressures)  and  spherical  geometries,  burning  rates  of  10  Ib/sec  translates  to  3.41  in. /sec,* 

7 Ib/sec  to  3.03  in. /sec,  and  4 Ib/sec  to  2.51  in. /sec.  The  comp=:tfcr  program  was  also  modified 
to  allow  for  radiant  heat  loss  from  the  burning  propellant  and  radiation  loss  from  the  gas  in 
the  chamber.  Summaries  of  the  output  from  these  various  techniques  are  contained  in  the 
tables  that  follow.  The  basic  assumptions,  in  addition  to  the  combustion  and  gas  flow  equations, 
are  contained  at  the  top  of  each  table.  We  feel  that  the  most  accurate  results  are  those  which 
allow  for  radial  burning  rates  (especially  3.40  in. /sec)  and  at  least  some  radiant  energy  loss. 

We  estimate  that  the  Category  V shield  has  an  effective  vent  area  of  7.785  sq  ft  when  vented. 

The  effective  vent  area  was  determined  by  combining  the  vent  areas  of  each  component  inside 
the  side  and  roof  panels  by  the  equation 

1*11 

= — + — 4-...  + — (166) 

Oeff  a,  Oj  On 

where 

Oeff  = effective  vent  area 

«i , Wj . . .ttn  = vent  area  of  each  componert 


*Timc  for  burtUng  vjrloui  quintitiet  (uniiiilng  H^hnicsl  geomeuici)  cin  be  calcuUled  from 

(1.  ^ ^ 

U>.  0.0602  Ib/in.v 


where  M it  Ihe  matt  of  the  propellant  in  pounds,  and  r is  the  burning  rate  in  inches/aecond.  For  3.41  in./sec  burning  rate, 
one  can  calculate  that  10  lb  burnt  In  1.00  second.  100  lb  In  2.  IS  seconds,  and  1000  lb  in  4.64  seconds- 
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8^  a:  4.848^ 


(165) 


1/3 

8 (164) 
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The  side  panel  consists  of 

• a layer  of  angle  irons  with  a vent  area  of  273  in.^ , 

r 

! • a layer  of  angle  irons  with  a vent  area  of  485  in.^ , 

t 

£ 

\ • three  perforated  plates  with  each  plate  having  a vent  area  of  1410  in.^ , 

• four  layers  of  aluminum  mesh  with  each  layer  having  a vent  area  of  3089  in.* . 
Thus,  from  Eq.  ( 1 66),  the  effective  vent  area  of  a side  panel  is 

273  485  1410  3089 

where  is  the  effective  vent  area  of  a side  panel. 

I 

The  roof  panel  consists  of 

• a layer  of  angle  irons  with  a vent  area  of  154  in.*, 

• a layer  of  angle  irons  with  a vent  area  of  273  in.* , 

j • three  perforated  plates  with  each  plate  having  a vent  area  of  795  in.* , 

i • four  layers  of  aluminum  mesh  with  each  layer  having  a vent  area  of  1 742  in.* . 

Thus,  from  Eq.  ( 1 66),  the  effective  vent  area  of  a roof  panel  is 


I ^ roof p.n,l  " 1 54 273''  795  ^1742 

1 

i >‘^oofp.n.l  = 61.6in.* 

I 

I 

i where  A y is  the  effective  vent  area  of  a roof  panel. 

f 

! Accounting  for  eight  side  panels  and  four  roof  panels,  the  total  effective  vent  area  of  the 

I Category  V shield  is 

! . (8X109.3)4(4X61.63)  _ 

i Ay  ^ 7.785  ft 

) 144 

I Figure  1 2 was  taken  from  the  numerical  results  contauied  in  Table  XIV  and  contains 

plots  of  peak  pressure  versus  vent  area  for  a burning  rate  of  150  Ib/sec,  no  radiation  loss,  and 
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FIGURE  12-  MIO  PROPELLANT  IN  VENTED  CHAMBER  WITH  NO  RADIATION  LOSS 


TABLE  XIV.  MtO  PROPELLANT  IN 
NO  RADIATION 


Auumptions:  Ambient  Pressure  - 14.7 

Iniiiul  Ratio  uf  Specific  Heats  » 1.4 

Ambient  Temperature  = 

Chamber  Volume  « 919.6 

No  radiation  loss 
Constant  mass/time  burning  rate 


Run 

No. 

AT  (sec) 

Effective 

Vent  Area 
(sq  ft) 

Burning 

Rate 

(Ibs/sec) 

1 

0.001 

-0- 

150 

2 

0.001 

-0- 

150 

3 

0.001 

-0- 

150 

4 

0.001 

0.006 

50 

S 

0.001 

0.006 

50 

6 

0.001 

0.006 

50 

7 

0.001 

0.006 

50 

8 

0.001 

0.006 

100 

9 

0.001 

0.006 

150 

10 

0.001 

0.006 

150 

11 

O.OOl 

0.006 

150 

12 

0.001 

0.006 

ISO 

13 

0.001 

0.006 

300 

14 

0.001 

0.18 

7 

15 

0.001 

0.18 

10 

16 

0.001 

1.2 

ISO 

17 

0.001 

1.2 

150 

18 

0.001 

1.2 

150 

19 

0.001 

150 

20 

0.001 

150 

21 

0.001 

ISO 

22 

0.001 

4.71 

10 

23 

0.001 

4.71 

10 

24 

0.001 

4.71 

10 

25 

0.001 

4.71 

SO 

26 

0.001 

4.71 

50 

27 

0.001 

4.71 

50 

28 

0.001 

4.71 

50 

29 

0.001 

4.71 

100 

30 

0.001 

4.71 

100 

31 

0.001 

4.71 

100 

32 

0.001 

4.71 

100 

33 

0.0005 

4.71 

124.4 

34 

0.001 

4.71 

124.4 

35 

0.005 

4.71 

124.4 

36 

0.01 

4.71 

124.4 

37 

0.001 

4.71 

124.4 

38 

0.001 

4.71 

124.4 
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LOSS 


CHAMBER  WITH 


psi 


cu  ft 


Quantity 

(lbs) 

Peak 

Pressure 

(psi) 

Remarks 

10 

86 

@ 0.067  sec 

50 

372 

@ 0.100  sec 

100 

535 

@ 0.334  sec 

10 

86 

@ 0.200  sec 

30 

235 

@ 0.600  sec 

50 

>368 

max.  not  reached 

100 

>368 

max.  not  reached 

30 

236 

@ 0.300  sec 

10 

86 

@ 0.067  sec 

30 

236 

@ 0.200  sec 

50 

371 

@ 0.334  sec 

100 

530 

@ 0.667  sec 

30 

237 

@ 0.100  sec 

10 

48 

max.  not  reached 

10 

67 

max.  not  reached 

10 

77 

@ 0.067  sec 

50 

173 

@ 0.255  sec 

100 

173 

@0.255  sec 

10 

69 

@ 0.067  sec 

50 

108 

@0.203  sec 

100 

108 

@ 0.203  sec 

10 

3 

@ 0.120  sec 

50 

3 

@0.120  sec 

100 

3 

@0.120  sec 

10 

24 

@0.136  sec 

30 

24 

@0.136  sec 

50 

24 

@0.136  sec 

100 

24 

@0.136  sec 

10 

44 

@ 0.100  sec 

30 

45 

@0.121  sec 

50 

45 

@0.121  sec 

100 

45 

@0.121  sec 

10 

51 

@ 0.080  sec 

10 

51 

@ 0.080  sec 

10 

49 

@ 0.080  sec 

10 

46 

@ 0.080  sec 

50 

54 

@0.1 17  sec 

100 

54 

@0.1 17  sec 

TABLE  XIV.  MIO  PROPELLANT  IS  VENTED  CHAMBER  WITH 
NO  RADIATION  LOSS  ICont'il) 


10-,  50-,  and  100-lb  quantities.  Figure  13  was  taken  from  the  numerical  results  contained  in 
Table  XV  and  contains  a plot  of  peak  pressure  versus  vent  area  for  a radial  burning  rate  of 
3.41  in. /sec  and  radiation  loss  from  the  burning  surface  and  surrounding  gas.  Figure  13  is 
for  1 0-lb  quantities  and  probably  predicts  peak  overpressure  more  accurately  than  Figure  1 2. 
Note  that  in  the  tables  and  figures,  pressure  is  in  units  of  psig;  that  is,  pressure  above  atmos- 
pheric pressure. 

TABLE  .XV.  M 10  PROPELLANT  IN  VENTED  CHAMBER 
WITH  RADIATION  LOSS 


i 

! 


I 


Assumptions:  Ambient  Pressure  = 14.7  psi 

Initial  Ratio  of  Specific  Heats  1 .4 

Ambient  Temperature  - 77°F 

Chamber  Volume  = 919.6  cu  ft 

Radiation  from  produced  gas  and  surrounding  gas 
Constant  Radid  Burning  Rate 
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FIGURE  13.  MIO  PROPELLANT  IN  VENTED  CHAMBER  WITH  RADIATION  LOSS 


D.  Comparison  with  Experimental  Results 

Several  tests  have  been  conducted  at  NSTL  where  MIO  propellant  has  been  Ignited  inside 
a Category  V suppressive  structure.  Although  a written  report  has  not  yet  been  published 
Dr.  Gary  McKown  of  NSTL  has  revealed  to  us  in  a telephone  conversation*  that  uncon- 
I'ined  MIO  propellant  is  u slow  burner.  During  the  Category  V tests,  100-,  2S0-,  and  600-lb 
quantities  of  MIO  were  burned  inside  the  structure.  The  propellant  was  in  cardboard  boxes 
no  more  than  a foot  tall  and  was  ignited  from  the  top.  The  material  burned  like  a “slow 
candle,”  with  a burning  time  of  approximately  two  or  three  seconds,  depending  on  the 
weight  of  the  propellant.  It  also  appeared  that  the  burning  rate  varied  linearly  with  the 
mass  of  propellant. 

During  the  tests  there  were  large  fireballs,  no  detonations,  and  very  low  quasi-static 
pressures  inside  the  structure.  From  the  films  of  the  event,  it  appeared  as  if  either  the  pro- 
pellant or  partially  combusted  products  of  the  reaction  had  gotten  outside  the  suppressive 
structure  and  ignited,  thus  contributing  to  the  enormous  fireball.  Also,  during  the  tests, 
ail  of  the  aluminum  screening  inside  the  structure  and  portions  of  other  aluminum  components 
were  burned  out  or  deteriorated,  possibly  due  to  the  intense  heat  of  the  reaction.  For  one 
experiment,  the  top  of  the  structure  was  covered  with  a piece  of  steel  to  control  the  vertical 
flow  of  the  fireball.  This  steel  roof  had  a hole  burned  through  it  by  the  end  of  the  tests, 

The  maximum  pressure  recorded  during  the  tests  was  approximately  0.2  psig.  Most  of  the 
time,  however,  pressure  was  so  small  that  it  was  difficult  to  record. 

As  noted  above,  we  suggested  that  Figure  13,  with  a radial  burning  rate  of  3.41  in. /sec 
which  reduces  to  100  lb  burned  in  2.15  seconds,  250  lb  burned  in  2.92  seconds,  and  600  lb 
burned  in  3.92  seconds,  was  the  best  figure  to  use  to  determine  peak  quasi-static  pressure. 

From  Table  XV  and  Figure  13,  for  an  effective  vent  area  of  4.71  square  feet  and  either  10, 

50,  or  100  lb  of  MIO  propellant,  one  can  see  that  the  peak  pressure  is  1 .4  psig.  The  initial 
Category  V suppressive  structure  has  an  estimated  effective  vent  area  of  7.785  square  feet. 
Since  the  interior  aluminum  baffles  deteriorated  after  a few  tests,  the  effective  vent  area 
increased  16.5%  to  9.069  ft^ . This  new  effective  vent  area  was  determined  by  eliminating 
the  last  term  in  the  equations  used  to  calculate  the  effective  vent  areas  of  the  side  panels 
and  roof  panels.  When  one  considers  the  effect  of  this  additional  vent  area  and  the  effect 
of  energy  losses  due  to  burning  of  propellant  outside  the  chamber,  one  would  expect  that 
the  curve  in  Figure  13  is  an  upper  limit  conservative  calculation.  If  these  other  terms,  that 
is,  additional  vent  area  and  radiant  heat  loss,  were  introduced  into  the  computer  program, 
lower  and  more  accurate  predictions  would  undoubtedly  have  been  produced. 


VI.  DISCUSSION 

Suppressive  shielding  was  developed  to  house  hazardous  operations  at  explosives 
processing  facilities.  Should  an  accidental  explosion  occur,  the  shield  is  designed  to  contain 
or  sharply  reduce  the  blast  overpressures, thermal  effects,  and  fragment  hazards  produced. 
Nearby  operations  are  thus  shielded  from  the  detrimental  effects  of  the  accident. 

To  properly  design  a shield  to  attenuate  or  completely  contain  the  products  of  the 
explosion  requires  a knowledge  of  several  complicated  phenomena.  These  include: 
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1 . overpressures  inside  and  outside  a vented  enclosure  produceu  by  the  detonation  of 
high  explosive  charges  or  by  the  burning  of  propellants, 

2.  fragments  generated  by  the  accident, 

3.  thermal  environment  produced  by  the  accident, 

4.  transient  loading  on  the  structural  components  of  the  shield  produced  by  the 
overpressures,  fragments,  and  thermal  environment. 

5.  elastic-plastic  response  of  the  shield  to  the  transient  loading, 

6.  resistance  of  the  shield  to  fragment  penetration,  and 

7.  attenuation  of  thermal  effects  by  the  shield. 

Other  aspects  of  the  design  include  problems  of  entry  and  exit  (for  per.  onnel,  explosives, 
and  utilities),  cleanup,  weatherproofing,  lighting,  air  conditioning,  and  economy. 

All  of  these  problems  have  been  addressed  during  the  suppressive  structures  program, 
but  not  all  of  them  have  been  satisfactorily  resolved.  Prediction  of  the  thermal  environment 
'ind  its  attenuation  by  the  shield  has  received  the  least  emphasis.  Most  other  aspects  of  the 
‘.esign  have  been  investigated  in  some  detail,  and  procedures  suitable  for  sh  eld  design  have 
oeen  developed;  however,  airiest  all  aspects  of  the  design  can  benefit  from  additional  research. 

The  work  documented  in  this  report  addresses  principally  items  1 , 4,  and  5.  This  work 
has  included  further  developments  in  the  application  of  energy  methods  for  the  analysis  of 
structural  response  to  blast  loading,  comparisons  between  predictions  of  structural  response 
and  test  results,  and  predictions  of  pressure-time  histories  inside  vented  enclosures  produced 
by  burning  propellant.  The  procedures  developed  for  predicting  pressures  inside  vented 
enclosures  produced  by  a burning  propellant  are  applicable  to  the  design  of  Category  V 
shields.  Formulas  and  graphical  solutions  for  structural  analysis  are  applicable  to  the  limit 
design  of  all  shield  groups. 

Comparisons  between  analysis  and  experiment,  which  are  described  in  this  report  and 
in  earlier  work*  * * have  shown  that  approximate  energy  methods  pro-'de  good  estimates 

for  the  elastic  and  elastic-plastic  behavior  of  simple  structural  elements  subjected  to  blast 
loading.  This  report  contains  a comprehensive  summary  of  solutions  developed  during  the 
suppiessive  structures  program  for  the  analysis  of  a large  number  of  different  structural 
elements.  The  comparisons  also  show  that  solutions  for  simple  structural  elements  are  not 
always  suitable  for  the  analysis  of  complex  structures.  Coupling,  which  can  occur  between 
structural  components,  can  greatly  attenuate  the  response  of  each  component  treated 
separately.  However,  even  in  this  case  we  have  shown  that  approximate  energy  solutions 
can  give  good  results,  but  with  some  increase  in  complexity.  Energy  methods  also  lend  then.- 
selves  to  graphical  solutions  which  greatly  simplify  the  complexity  and  labor  involved  in  the 
solution  process. 
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Although  substantial  progress  has  been  made  in  the  development  of  analytical  procedures 
to  support  the  design  of  suppressive  shields  additional  work  is  needed  in  several  areas.  These 
include: 

1.  definition  of  the  thermal  environment  produced  by  the  detonation  of  a high 
explosive  or  by  the  burning  of  propellant, 

2.  methods  to  predict  the  thermal  environment  outside  of  a suppressive  shield, 

i.  comparisons  between  measured  and  predicted  pressures  produced  by  burning 

propellants  in  vented  enclosures, 

4.  failure  criteria  for  dynamically  loaded  commercially  available  structural  components. 
Commercially  available  components  have  been  designed  for  optimum  efficiency 
under  static  loading.  A better  understanding  of  their  behavior  when  loaded  by 
short  duration  high  intensity  loads  is  needed.  This  is  particularly  true  for  the 
failure  of  wide-flange  beams  in  shear. 
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